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Abstract 



We perform a complete analysis of the supersymmetric particle spectrum in the Min- 
imal Supergravity (mSUGRA) model where the soft SUSY breaking scalar masses, 
gaugino masses and trilinear couplings are unified at the GUT scale, so that the elec- 
troweak symmetry is broken radiatively We show that the present constraints on the 
Higgs boson and superparticle masses from collider searches and precision measure- 
ments still allow for large regions of the mSUGRA parameter space where charginos, 
neutralinos, sleptons and top squarks as well as the heavier Higgs particles, are light 
enough to be produced at the next generation of e + e~ linear colliders with center of 
mass energy around y^ ~ 800 GeV, with sizeable cross sections. An important part 
of this parameter space remains even when we require that the density of the lightest 
neutralinos left over from the Big Bang, which we calculate using standard assump- 
tions, falls in the range favored by current determinations of the Dark Matter density 
in the Universe. Already at a cm. energy of 500 GeV, SUSY particles can be accessible 
in some parameter range, and if the energy is increased to yfs ~ 1.2 TeV, the e + e~ 
collider will have a reach for high precision studies of SUSY particles in a range that 
is comparable to the discovery range of the LHC. 



1. Introduction 

Supersymmetric theories (SUSY) [|TJ are the best motivated extensions of the Standard Model 
(SM) of the electroweak and strong interactions. They provide an elegant way to stabilize the 
huge hierarchy between the Grand Unification (GUT) or Planck scale and the Fermi scale, 
providing a natural framework to cancel the quadratic divergences of the radiative corrections 
to the Higgs boson mass. The most economical low-energy supersymmetric extension of 
the SM, the Minimal Supersymmetric Standard Model (MSSM) [0, allows for a consistent 
unification of the three coupling constants of the SM gauge group |§. In addition, it can 
provide a natural solution of the Dark Matter problem 0], since it predicts the existence 
of an electrically neutral, weakly interacting, massive and absolutely stable particle; for 
large regions of parameter space the thermal relic density of this particle agrees with the 
Dark Matter density derived from cosmological arguments. The search for Supersymmetric 
particles and for the required extended Higgs spectrum is one of the main motivations for 
building high-energy colliders. 

In the MSSM one assumes the minimal gauge group, i.e. the SM SU(3)c x SU(2) L x U(l)y 
group; the minimal particle content, i.e. three generations of fermions [without right-handed 
neutrinos] and their spin-zero partners as well as two Higgs doublet superfields to break the 
electroweak symmetry; and R-parity conservation || which makes the lightest SUSY parti- 
cle, assumed to be the lightest neutralino x?, absolutely stable. In order to explicitly break 
supersymmetry [as required by experiment] while preventing the reappearance of quadratic 
divergences, a collection of soft terms is added to the Lagrangian ||: mass terms for the 
gauginos, mass terms for the scalar fermions, mass and bilinear terms for the Higgs bosons, 
and trilinear couplings between sfermions and Higgs bosons. In the general case, that is if one 
allows for intergenerational mixing and complex phases, the soft SUSY breaking terms will 
introduce a huge number (105) of unknown parameters |7j], in addition to the 19 parameters 
of the SM. 

This feature makes any phenomenological analysis in the general MSSM a daunting task, 
if possible at all. In addition, almost all "generic" sets of these parameters are excluded by 
severe phenomenological constraints, on flavor changing neutral currents (FCNC), additional 
CP-violation, color and charge breaking minima, etc. Almost all FCNC problems are solved 
at once if the MSSM parameters obey a set of universal boundary conditions at the GUT 
scale. We will take these parameters to be real, which solves all problems with CP violationF]. 
The underlying assumption is that SUSY-breaking occurs in a hidden sector which commu- 
nicates with the visible sector only through gravitational-strength interactions, as specified 

1 If soft breaking parameters are universal at the GUT scale, they are allowed to have large CP-violating 
phases only in certain very narrow regions of parameter space where large cancellations occur between 
various contributions to electric dipole moments of the electron and neutron || . With all high-scale soft 
breaking parameters being real, the model predicts very small deviations from the SM in B— meson mixing 
and CP-violating decays || now being explored at the B— factories. However, we will see later that the 
model allows for significant new contributions to b — > sj decays, and related decay modes. 

1 



by supergravity [10]. Universal soft breaking terms then emerge if these supergravity inter- 
actions are "flavor-blind" [like ordinary gravitational interactions]. This is assumed to be 
the case in the constrained MSSM or minimal Supergravity (mSUGRA) model ]TT|. In this 



model the entire spectrum of superparticles and Higgs bosons is determined by the values of 
five free parameters. Since universal boundary conditions imply that the electroweak sym- 
metry is broken radiatively, which imposes one constraint on the input parameters, one is 
left with only four continuous free parameters and a discrete one. This makes comprehensive 
scans of parameter space feasible. 

Although other viable SUSY models exist, the mSUGRA model has become the most 
frequently used benchmark scenario for supersymmetry, and has been widely used to analyze 
the expected SUSY particle spectrum and the properties of SUSY particles, and to compare 
the predictions with available and/or expected data from collider experiments. Several global 
or partial analyses of the present theoretical and experimental constraints on the mSUGRA 



model have been performed in the literature; see for instance Refs. []12| , |13[ [14 . 

In this paper, we perform an independent analysis of the SUSY particle spectrum in the 
mSUGRA model, taking into account theoretical constraints^ and all available experimental 



information fL5| : searches for the MSSM Higgs bosons and SUSY particles at the LEP and 
Tevatron colliders, electroweak precision measurements, the radiative b — ► sj decay, etc. 
Special attention is paid to the implications of the measurement of the anomalous magnetic 
moment of the muon recently performed at Brookhaven |16| , and to the ~ la evidence for 
a SM-like Higgs boson with a mass Mmggs ~ 115.6 GeV seen by the LEP collaborations 
H] . We also discuss the implication of requiring thermal relic neutralinos to form the Dark 
Matter in the Universe. 

We show that in a large part of the mSUGRA parameter space at least one of these 
independent pieces of evidence for physics beyond the SM, from the Dark Matter density, the 
(g M — 2) measurement and the LEP2 Higgs boson-like excess, can be explained in mSUGRA. 
On the other hand, only a small area of parameter space allows for these three constraints 
to be fulfilled simultaneously. If all these indications survive further scrutiny, the parameter 
space of the model would thus already be tightly constrained. However, one has to keep 
in mind that the statistical significance for the LEP Higgs signal and the (g^ — 2) anomaly 
are at present still quite weakP], while the calculation of the thermal relic density relies on 
additional assumptions that cannot be tested in collider experiments. 

We then discuss prospects for producing SUSY particles and the heavier Higgs bosons 

2 In order to further limit the parameter space, one could require that the SUGRA model is not fine-tuned 
and the SUSY breaking scale should not be too high, a constraint which can be particularly restrictive since 
sparticles with masses beyond ~ 1 TeV would be excluded. However, the degree of fine-tuning which can be 
considered acceptable is largely a matter of taste, so for the most part we disregard this issue in our analysis. 

3 Recent estimates of the uncertainties in the hadronic contributions to (g^ — 2) might slightly push the 
theoretical prediction 118] towards the SM value and thus decrease the significance of the discrepancy, see 
for instance Ref. flil. 



of the MSSM at future high-energy e + e linear collider^] with center of mass energies \fs 
around 800 GeV as expected, for instance, at the TESLA machine ||24j|. We show that in 



large areas of the mSUGRA parameter space the production rates for the lightest charginos 
and neutralinos, as well as for sleptons, top squarks and the heavier Higgs bosons are large 
enough for these particles to be discovered, given the very large integrated luminosities, 
£ ~ 500 fb _1 , expected at this collider. Even at lower energies, -Js ~ 500 GeV, charginos, 
neutralinos and tau sleptons can be produced in some parameter range. If the energy is 
raised to */s ~ 1.2 TeV ||25| , the e + e~ collider will have a reach for probing the SUSY 



particle spectrum and the heavy Higgs bosons which is comparable to the reach of the LHC. 
The remainder of this paper is organized as follows. In the next section we briefly 
summarize the main features of the mSUGRA model and the way it is implemented in 
our analysis. In section 3 the experimental and cosmological constraints on the mSUGRA 
parameter space are discussed. In section 4 we analyze the production of SUSY particles 
and MSSM Higgs bosons at high energy e + e~ colliders. Conclusions are given in section 
5. For completeness, expressions for the cross sections of all discussed particle production 
channels in e + e~ collisions are collected in the Appendix. 

2. The Physical Set-Up 

2.1 The mSUGRA model 

We will perform our analysis in the constrained MSSM or minimal Supergravity model, 
where the MSSM soft breaking parameters obey a set of universal boundary conditions at 
the GUT scale, so that the electroweak symmetry is broken radiatively. For completeness and 
to fix the notation, let us list these unification and universality hypotheses and summarize 
the main features of the radiative electroweak symmetry breaking (EWSB) mechanism. 

Besides the unification of the gauge coupling constants, which is verified given the ex- 
perimental results from LEP1 || and which can be viewed as fixing the Grand Unification 
scale Mgut ~ 2 • 10 16 GeV, the unification conditions are as follows: 

- Unification of the gaugino [bino, wino and gluino] masses: 

Mi(M G ut) = M 2 (M GVT ) = M 3 (M GXJT ) = m 1/2 . (1) 

- Universal scalar [sfermion and Higgs boson] masses [i is a generation index] 

M Qi (M GVT ) = i\4 k (M GUT ) = M d - k (M GUT ) = M Zl (M GUT ) = M^(M G ut) 

= M Hu (M GVT ) = M Hd (M GVT ) = m . (2) 

4 During the final stage of the present work, for which preliminary results have been presented in Ref. p0| , 
the paper "Proposed Post-LEP Benchmarks for Supersymmetry" |l3| , which discusses some of the issues 
considered here appeared on the web. A brief discussion of the prospects at future e + e~ colliders has 
also been given in Ref. pl| . Prospects for the Tevatron Run II and for the LHC have been discussed in 
Refs. |^, [23], respectively; of course, these earlier studies used slightly weaker experimental constraints, as 
was appropriate at the time of their writing. 



- Universal trilinear couplings: 

A$(M GTJT ) = ^'(Mqut) = A?'(M GUT ) = 4)V (3) 

Besides the three parameters mi/ 2 , wi and v4 , the supersymmetric sector is described at 
the GUT scale by the bilinear coupling B and the supersymmetric Higgs(ino) mass parameter 
ji. However, one has to require that EWSB takes place. This results in two necessary 
minimization conditions of the two Higgs doublet scalar potential which, at the tree-level, 
has the form |26| [to have a more precise description, one-loop corrections to the scalar 
potential have to be included, as will be discussed later]: 

^Higgs = mlHlH d + rn 2 2 HlH u + ml{H u -H d + h.c.) 

+ 9 -^(HlH d -HtH u r + 9 l(HlH u )(HlH d ), (4) 

where we have used the usual short-hand notation: 

— 2 2 2 — 2 2 , 2 — 2 n /c-\ 

m l = m H d + V » m 2 = m H u + V , m 3 = Bfi. (5) 

The SU(2) invariant product of two doublets is defined as 0i-02 = 0i02 — 0i02> where the su- 
perscripts are SU(2) indices. The two minimization equations dVm g g S /dH d = dVm g g S /dH® = 
can be solved for /i 2 and B\i: 

1 



tf 



Bfi = - sin 2/3 



tan 2(3 (m H tan (3 — m H cot/5) — M z 
m H u + m H, + 2 /i : 



(6) 



Here, M\ = (gf + g%) ■ (t> 2 + vj)/4, and tan/3 = v u /v d is defined in terms of the vacuum 
expectation values of the two neutral Higgs fields. Consistent EWSB is only possible if 
eq. (||) gives a positive value of /i 2 . The sign of /i is not determined. Therefore, in this model 
one is left with only four continuous free parameters, and an unknown sign: 

tan/3 , mi/2 , ra , A , sign(ju). (7) 

All the soft SUSY breaking parameters at the weak scale are then obtained through Renor- 
malization Group Equations (RGE) j27|. 



The number of parameters could be further reduced by introducing an additional con- 
straint which is based on the assumption that the b and r Yukawa couplings unify at the 
GUT scale, as predicted in minimal SU(5). This restricts tan/3 to two narrow ranges around 
tan/3 ~ 1.5 and tan/3 ~ mt/rrib [01 • The low tan/3 solution is ruled out since it leads to a 
too light an h boson, in conflict with searches at LEP2. However, Yukawa unification is not 
particularly natural in the context of superstring theories, and minimal SU(5) predictions are 
known to fail badly for the lighter generations. We therefore treat all three third generation 
Yukawa couplings as independent free parameters. 



2.2 Calculation of the SUSY particle spectrum 

In this section, we briefly discuss our procedure for calculating the SUSY particle spectrum 
in the constrained MSSM with universal boundary conditions at the GUT scale, as well 
as related issues which are relevant to our study. All results are based on the numerical 



FORTRAN code SuSpect version 2.0 p9| , to which we refer for a more detailed description. 
The algorithm essentially includes: 

- Renormalization group evolution (RGE) of parameters between the low energy scale 
[Mz and/or the electroweak symmetry breaking scale] and the GUT scale. 

- Consistent implementation of radiative electroweak symmetry breaking (EWSB). Loop 
corrections to the effective potential are included using the tadpole method. 

- Calculation of the physical (pole) masses of the Higgs bosons, scalar quarks and leptons 
as well as gluinos, charginos and neutralinos. 

In more detail we proceed as follows. We first chose the low-energy input values of the SM 
parameters. The gauge couplings constants are defined in the MS scheme at the scale Mz 



_ Aira^jMz) 2 _ Ana^(M z ) , _ m , 

9\ - =2 ' 92 - =2 > 93 - ^ a s \ M Z, 



Their values have been obtained from precision measurements at LEP and Tevatron fT5|| : 
a^(M z ) = 1/127.938 , af^(M z ) = 0.1192 , s 2 w = 0.23117. (9) 



The pole masses of the heavy SM fermions are [15] 



M t = 174.3 GeV , M b = 4.62 GeV , M T = 1.778 GeV . (10) 



From the pole 6-quark mass, one then obtains the DR mass, m^rrTf,) ~ 4.23 GeV which 
is then evolved, using two-loop 0(a 2 s ) RGE, to obtain the running mass at the scale M z , 
mb(Mz) — 2.92 GeV. Since the two-loop corrections to the difference between pole and DR 
top and bottom quark masses are not yet known, we include, instead, the analogous two- 
loop corrections in the MS scheme, which should be close to the DR ones. The difference 
should not be important in view of the experimental errors in the determination of the two 



masses |Lj], M t = 174.3 ± 5.1 GeV and m b (m b ) ~ 4.24 ± 0.11 GeV. 



Next, the DR-scheme values of the gauge and Yukawa couplings are extracted from these 



inputs [|30 1 . The latter are defined by [v = 174.1 GeV]: 



MMt ) = ^ , MMz) = ^ma , MMz) = ^m. . (11) 

v sin p v cos p v cos p 



All couplings are then evolved up to the GUT scale using two-loop RGEs |3Q, 31 . Here 



heavy (super)particles are taken to contribute to the RGE only at scales larger than their 



mass, i.e. multiple thresholds are included in the running of the coupling constants near 
the weak scale. The GUT scale Mqut — 2 • 10 16 GeV is defined to be the scale at which 



9i — 92 • v3/5. We do not enforce #2 = gz at the GUT scale and assume that the small 
discrepancy [of the order of a few percent] is accounted for by unknown GUT-scale threshold 
corrections [ |32|| . 

In our numerical analyses we fix the MSSM parameters tan/? [given at scale Mz] as well 
as Ao and the sign of /z, and then perform a systematic scan over the high energy mSUGRA 
inputs mo and mm. Given these boundary conditions, all the soft SUSY breaking parameters 
and couplings are evolved down to the electroweak scale. Our default choice for this scale is 
the geometric mean of the two top squark masses, Mewsb = ij^i-^h) 1 ^ 2 •> w hich minimizes 
the scale-dependence of the one-loop scalar effective potential [[33]]. Since tan/5 is defined at 



scale Mz, the vevs have to be evolved down from M EWSB to M z p3|j . 

One-loop radiative corrections to the Higgs potential play a major role in determining 
the values of the parameters |/i| and B in terms of the soft SUSY breaking masses of the 
two Higgs doublet fields. We treat these corrections using the tadpole method. This means 
that we can still use eq. @ to determine /x 2 (M E wsb); one simply has to add one-loop 



tadpole corrections to m H and mjj @, pj|. We include the dominant third generation 



fermion and sfermion loops, as well as subdominant contributions from sfermions of the first 
two generations, gauge bosons, Higgs bosons, charginos and neutralinos, with the running 
parameters evaluated at M EW sb- As far as the determination of /z 2 and B\i is concerned, this 
is equivalent to computing the full one-loop effective potential at scale M EWSB . Since \\i\ 
and B affect masses of some (s)particles appearing in these corrections, this procedure has 
to be iterated until stability is reached and a consistent value of \x is obtained; usually this 
requires only three or four iterations for an accuracy of O(10~ 4 ), if one starts from the values 
of \\x\ and B as determined from minimization of the RG-improved tree-level potential at 
scale M EW sb- 

At this stage, we check whether the complete scalar potential has charge and/or color 
breaking (CCB) minima, which can be lower than the electroweak minimum. These can e.g. 
appear in the top squark sector]] for large values of the trilinear coupling A t . In order to 
avoid them, we impose the (simplest) condition [57]: 



CCB: A 2 t <3(ml+ml+^ + m 2 Hu ). (12) 

Of course, we also reject all points in the parameter space which lead to tachyonic Higgs 
boson or sfermion masses^: 

NoTachyon: M 2 A > , Ml > , m) > 0. (13) 



5 CCB minima involving first and second generation sfermions are usually separated from the desired 
EWSB minimum by high potential barriers, so that the EWSB minimum is still stable on cosmological time 
scales p6| . 

6 Later on, we will be more restrictive and discard the situations where SUSY particles have masses which 
are lower than the mass of the neutralino Xi which will be assumed to be the lightest SUSY particle 



The EWSB mechanism is assumed to be consistent when all these conditions are satisfied. 

We then calculate all physical particle masses. The procedure is iterated at least twice 
until stability is reached, in order to take into account: (i) Realistic (multi-scale) particle 
thresholds in the RG evolution of the dimensionless couplings via step functions in the /3 
functions for each particle threshold, (ii) Radiative corrections to SUSY particle masses, 
using the expressions given in Ref. j3"5f , where the renormalization scale is set to M EW sb- 

We first evaluate the SUSY-radiative corrections to the heavy fermion masses, rh t , fh^ 
and rh T , following Ref. |[I|.[] This includes SUSY-QCD corrections for the t,b quarks 
[from squark-gluino loops] and the dominant electroweak corrections for the b and r masses 
[chargino-sfermion loops which are enhanced by terms oc /xtan/3]. As suggested in Ref. fl38| . 
we use the "MSSM" quark masses [essentially the Yukawa coupling times vev] in the squark 
mass matrices. Our iteration then resums all SUSY-QCD corrections to the quark masses 
of order (a s tan/3) n . This is important at large tan/3, where these corrections can be quite 
sizable. The various sectors of the MSSM are then treated as follows: 

- In the sfermion sector, the soft scalar masses as well as the trilinear couplings for the 
third generation are obtained using one-loop RGE, and are frozen at the scale M EW sb- In 
the third generation sfermion sector [t, b, f], mixing between "left" and "right" current eigen- 
states is included, where the radiatively corrected running fermion masses at scale M EWSB 
are employed in the sfermion mass matrices. The radiative corrections to the sfermion masses 
are included according to Ref. P5| , i.e. only the QCD corrections for the superpartners of 
light quarks [including the bottom squark] plus the leading electroweak corrections to the 
top squarks; the small electroweak radiative corrections to the slepton masses have been 
neglected. 

- In the gaugino sector, the SUSY breaking gaugino masses are obtained using the two- 
loop RGEs and are also frozen at M E wsb- The mass matrices for charginos and neutralinos 
are diagonalized using analytical formulae |3J|. The one-loop QCD radiative corrections to 



the gluino mass are incorporated [p0|| ) while in the case of charginos and neutralinos the 
radiative corrections |4(| are included in the gaugino and higgsino limits, which is a very 



good approximation according to Ref. [35 



- In the Higgs sector, the running mass of the pseudoscalar Higgs boson is obtained from 
the soft SUSY breaking Higgs masses [again frozen at M EWSB ] and the full one-loop tadpole 
corrections |3J|. This mass is then used as input, together with M t , tan/3 and some MSSM 
parameters [A t , A b , \i and the soft SUSY breaking third generation squark masses], to obtain 
the pole masses of the pseudoscalar Higgs boson A, the two CP-even h and H bosons and 



the charged H Higgs particles. This last step is similar to the program HDECAY |41|], which 



calculates the Higgs spectrum and decay widths in the MSSM. The complete radiative cor- 



7 In our procedure some of the leading logarithmic terms are already included in the RG evolution of the 
Yukawa couplings via the step functions. Therefore, care has to be taken to avoid double counting when 
extracting the relevant radiative corrections from the expressions given in Ref. [|35[. 



rections due to top/bottom quark and squark loops within the effective potential approach, 
leading NLO QCD corrections [through renormalization group improvement] and the full 
mixing in the stop and sbottom sectors are incorporated using the analytical expressions of 



Ref. [|2|. We have verified that the results obtained for the Higgs spectrum, in particular for 
the lightest h boson mass, are nearly the same as those obtained from the complete results 
of the Feynman diagrammatic approach implemented in the program FeynHiggs [33|] . 

Our results for some representative points of the mSUGRA parameter space have been 
carefully cross-checked against other existing codes. We obtain very good agreement, at the 
one percent level, with the program SOFTSUSY f|4| which has been released recently^. We 
also find rather good agreement for the SUSY particle masses computed by the program 



ISASUGRA [ 451 , once we chose the same configuration [soft SUSY breaking masses frozen 
at Mz, some radiative corrections to sparticle masses are not included, etc..]. The value 
of the lightest Higgs boson mass is in less good agreement, presumably due to the more 
sophisticated treatment of the Higgs potential in SuSpect; we will see later that a precise 
calculation of the h boson mass is an important ingredient of our analysis. 

3. Constraints on the mSUGRA parameter space 

3.1 Experimental Constraints 

i) Lower bounds on the SUSY particles masses 

A wide range of searches for SUSY particles has been performed at LEP2 and at the Tevatron, 
resulting in limits on the masses of these particles [[HJ. The pair production of the lightest 
chargino at LEP2, e + e~ — > XiXii would probably have been the cleanest SUSY process. 
In general it has the largest SUSY production cross section at e + e~ colliders, after the 
experimental cuts needed to suppress the backgrounds, and the information that it provides 
is one of the most important in the context of the mSUGRA model. The negative outcome 
of searches for charginos at LEP2, up to energies of \/s ~ 208 GeV, gives the approximate 
bound m k ± > 104 GeV [0.0 



In mSUGRA the gaugino masses are unified at the GUT scale, leading to the approximate 
relation M 2 — 1M\ ~ M 3 /3 at the weak scale. The bound on the lightest chargino mass thus 
translates into a lower bound on the LSP mass, m^o > 50 GeV [in the gaugino-like region; 

8 We thank Ben Allanach for his gracious help in performing this detailed comparison. 

9 This bound is not valid if |/x| <C M2, i.e. for a very higgsino-like chargino which is almost degenerate 
with the LSP leading to a small release of missing energy. We therefore exclude slightly too much in the 
"focus point" region, see below. The true bound is also reduced somewhat in scenarios with light sneutrino, 
since v exchange in the t— channel reduces the cross section, and since \i — » v + ^ ± decays can be difficult 
to detect if the X\ — v mass difference is small; this can happen for very small too in mSUGRA, but such 
scenarios are tightly constrained by slepton searches and SUSY loop effects. An accurate treatment of this 
bound is not important for the main topic to be investigated in this paper, the reach for future e + e~ colliders 
with energy much above the LEP range. 



in the higgsino-region, the bound is higher] and also on the gluino mass, m g ~ M 3 > 300 
GeV. In the case of the LSP, the bound can be improved by using searches for neutralino 
production at LEP2, e + e _ —>■ XiX2'i however, these neutralino searches are relevant only 
for low values of tan (3 < 2 which are already excluded by Higgs boson searches, as will be 
discussed later. In the case of gluinos, this indirect bound is similar to the one obtained 
from direct searches at the Tevatron, m g > 300 GeV, which is valid if rriq ~ m g ||47|| ; the 
direct search limits for rriq 3> m g are significantly weaker. 

The bound on M2 also translates into a bound on the masses of first and second generation 
squarks. The RG evolution of these masses [up to small contribution from the D-terms] gives 
the approximate relation m~ ~ m\ + Qm\, 2l which leads to rriq > 250 GeV, again of similar 



size as the bound from direct searches at the Tevatron |fJ7 |. For third generation squarks 
the RGE are more complicated and mixing between the eigenstates is important, due to the 
large values of the Yukawa couplings, so that the bounds from direct searches are relevant. 



We use the bound from LEP2 [46|, which is almost independent of the decay modes, and is 
applicable down to squark— Xi mass splittings of a few GeV; Tevatron search limits jPJ are 
stronger in some cases, but more dependent on details of squark decay, and disappear for 
mass splittings below ~ 30 GeV. 

By far the tightest slepton search limits also come from LEP2 [fffi| . Here the coefficients 
of the m\ j 2 term in the RG evolution of m~ are small so that [at least for the SU(2) singlet 
"right-handed" states] mj ~ rriQ. These LEP bounds are generally a few GeV below the 
kinematical limit, except for some small regions of the parameter space with small mass 
difference to the LSP.Q Since sneutrinos might decay invisibly into vx\i only indirect bounds 
can be placed on their masses. Limits from searches for charged SU(2) doublet sleptons, 
whose masses are related to sneutrino masses by SU(2), are rather model-independent. In 
mSUGRA additional indirect limits on sneutrino masses follow from searches for SU(2) 
singlet sleptons and charginos, via the resulting constraints on mo and mm- 

To summarize, in our analysis we impose the following bounds on sparticle masses: 

m-+ > 104 GeV , m?> 100 GeV with / = i u b u f* v 
m g > 300 GeV , m^ 2 > 260 GeV with q = u,d, 5, d. (14) 

ii) Constraints from the Higgs boson masses 

The search for SUSY Higgs bosons was the main motivation for extending the LEP2 energy 



up to y/s ~ 208 GeV. In the SM, a lower bound [g|, Hi 



M H o > 113.5 GeV at the 95%C.L. (15) 



10 Contrary to chargino pair production, the cross sections for sfermion pair production in e + e collisions 
is suppressed by a /3 3 factor near threshold so that it is rather tiny near the kinematical limit. 



has been set by investigating the Higgs-strahlung process, e + e~ — > ZH°. In the MSSM, this 
bound is valid for the lightest CP-even Higgs particle h if its coupling to the Z boson is 
SM-like, i.e. if g'zzh/d'zzH — sin 2 (/3 — a) ~ 1 where a is the mixing angle in the CP-even 
neutral Higgs sector. This occurs in the decoupling regime where M\ 3> M§. For much 
lower values of Ma and for large tan/3, one has Mh ~ Ma and sin 2 (/3 - a) < 1. In this 
case, the bound (|15|) applies to the mass of the heavy neutral CP-even Higgs boson H. 
However, tighter constraints on the parameter space can usually be derived from the search 
for the associated production of CP-even and CP-odd Higgs particles, e + e _ —* hA, which 
has a cross section proportional to cos 2 (/3 — a). At LEP2, a combined mass exclusion in the 
MSSM, 

M h ~ M A > 93.5 GeV (16) 

at the 95% Confidence Level, has been set in this case. This limit is lower than that from 
the Higgs-strahlung process, due to the less distinctive signal and the /3 3 suppression near 
threshold for spin-zero particle pair production. 

Deriving a precise bound on Mh for arbitrary values of Ma [not just for M\ 3> M§ or 
Ma — Mh] and hence for all values of sin 2 (/3 — a) is rather complicated, since one has to 
combine results from two different production channels, which have different kinematical 
behavior, cross sections, backgrounds, etc. In our analysis, we will use an interpolation 
formula for the excluded value of M^: 

M h > 93.5 + 15a; + 54.3a; 2 - 48.4a; 3 - 25.7a; 4 + 24.8a; 5 - 0.5 GeV, (17) 

with x = sin 2 (/3 — a). This formula fits the exclusion plot in the [Mh, sin 2 (/? — a)} plane 



given by the ALEPH collaboration^ in Ref. [50|. At this point we assume a very small 
theoretical error of 0.5 GeV in the calculation of Mh and Ma, and have lowered the bounds 
accordingly; this is the typical difference between our value of Mh and the one obtained in 
the full diagrammatic approach |i3" ]. The total theoretical error on the calculated Higgs 



boson masses is probably much larger, even if all input parameters were known exactly. The 
intrinsic uncertainty due to unknown higher-order effects is usually estimated to be about 
two or three GeV. In addition, the error on the top quark mass leads to an uncertainty of the 
predicted value of Mh of a few GeV.0 The impact of these uncertainties will be discussed 
later. 

We will also study the implications of the 2.1<r evidence for a SM-like Higgs boson with 
a mass Mu = 115.6 GeV seenQ by the LEP collaborations JT7]]. In view of the theoretical 



11 The ALEPH search was performed in the energy range up to y/s = 202 GeV leading to the bounds 
M h > 107.7 GeV in the decoupling limit and M h , M A > 91.5 GeV for sin 2 (/3 - a) - 0. We have extended 
these end points to the values Mh ~ 113.5 and 93.5 GeV obtained at y/s ~ 208 GeV, while leaving the 
general form of the exclusion contour unchanged. 

12 The finite gluino contributions to the two-loop radiative corrections J43|, which we did not take into 
account, might also change Mh by one or two GeV. 

13 This used to be interpreted as ~ 2.9er evidence for a Higgs boson with mass Mh — 115.0 GeV. However, 
a recent re-evaluation led to a lower significance, and correspondingly slightly higher preferred Higgs mass. 
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and experimental uncertainties, we interpret this result as favoring the range 

113 GeV < M h < 117 GeV and sin 2 (/5 - a) > 0.95. (18) 

The second requirement ensures that the cross section for Higgs production, a(e + e~ — > 
hZ), is similar to that of the SM [the excess is assumed to come from the Higgs-strahlung 
pro cess]. 

iii) Constraints from electroweak precision observables 

Loops of Higgs and SUSY particles can contribute to electroweak observables which have 
been precisely measured at LEP, SLC and the Tevatron. In particular, the radiative correc- 
tions to the self-energies of the W and Z bosons, Hww an d H-zz, might be sizeable if there 
is a large mass splitting between some particles belonging to the same SU(2) doublet; this 
can generate a contribution which grows as the mass squared of the heaviest particle. The 
dominant contributions to the electroweak observables, in particular the W boson mass and 



the effective mixing angle sfy, enter via a deviation of the p parameter |51| from unity, which 
measures the relative strength of the neutral to charged current processes at zero momentum 
transfer, i.e. the breaking of the global custodial SU(2) symmetry: 



< IV 



c w M w 2 2c w s w 



Ap , As' w ~ - 7 2 A P 



9(V 2 - s 2 ) F ' w ~ r 2 - s 2 

\-l . a rr rrW/n/r? rr fr\\/i\/r2 



p=(l- Ap)' 1 ; Ap = U zz (0)/M z - n ww (0)/M* r . (19) 



Let us briefly summarize the possible contributions in the MSSM [^, |52|, [53|. A close 



inspection of the contributions of the Higgs and chargino/neutralino sectors shows that 
they are very small, Ap < 10~ 4 . In the former case, the contributions are logarithmic, 
~ alog(Mh/Mz), and are similar to those of the SM Higgs boson [they are identical in the 
decoupling limit] . In the latter case they are small because the only terms in the chargino and 
neutralino mass matrices which could break the custodial SU(2) symmetry are proportional 
to My/. Since first and second generation sfermions are almost degenerate in mass, they 
also give very small contributions to Ap.0 Only the third generation sfermion sector can 
generate sizable corrections to the p parameter, because left-right mixing and [in case of the 
stop] the supersymmetric contribution oc rn\ leads to a potentially large splitting between 



the sfermion masses. In particular, the contributions of the (£, b) and (r, v) iso-doublets [53 



14 The "discovered" Higgs boson could also be the heavy CP-even H boson produced in the strahlung 
process with cos 2 (/3 — a) close to unity, if the h and A bosons have masses Mh ~ Ma > 93.5 GeV. However, 
in mSUGRA only a very small set of input parameters can give this configuration. In fact, we did not find 
any parameter set where this possibility is realized and all other constraints are satisfied. 

15 It has recently been pointed out 154] that light sleptons and light charginos could slightly improve the 
quality of global fits to electroweak precision data. However, such configurations are excluded in mSUGRA, 
chiefly by the Higgs boson search limits. 
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might become unacceptably large. They are given by: 

3Crp 



Ap(t,b) 



Ap{f, v) 



2 2 el 2 2 \ 2 2 r/ 2 2 \ 

8tt 2 a/2 L ° tSt ^ m ^ ,m ^ ~ ° bSb ^ m b^ m b 2 ) 
+ 4 [ c fe f( m h m t) + s l f( m h m i 2 )} + St [cl f{m\, m~ bi ) + s 2 b /(m|, m|)] 

Grp 



8<K 2 y/2 



-c T s T f(m fl ,m~ 2 ) + c T f(m fi ,m~) + s T f(m~ 2 ,m^) 

f(x,y) = x + y-2^-\og-, (20) 

x-y y 



where Si and q are the sine and cosine of the mixing angles. In the stop sector large values 
of the trilinear coupling A t can be dangerous, since they lead to values m^ <C m^ , m~ b . In 
the stau/sbottom sectors large values of the parameters tan (5 and \i lead to sizable splitting 
between the b, f mass eigenstates. 

We have required the contributions of the third generation sfermions to stay below the 
acceptable level of 

Ap(f) < 2.2 ■ 1(T 3 , (21) 



which approximately corresponds to a la deviation from the SM expectation ||55|| . Two-loop 
QCD corrections to 5p(t, b) have been calculated |)6] and can be rather large, increasing the 



contribution by up to 30%. In our analysis we include the leading components of these 
corrections, i.e. the full corrections due to gluon exchange and the correction due to gluino 
exchange in the heavy gluino limit [which should be a good approximation in our case]. 
However, in our numerical analyses we find that this constraint is usually superseded by the 
Higgs boson search constraint. 

iv) The b — > 57 constraint 

Another observable where SUSY particle contributions might be large is the radiative flavor 



changing decay b — > 57 ]57||. In the SM this decay is mediated by loops containing charge 
2/3 quarks and iy-bosons. In SUSY theories additional contributions come from loops 
involving charginos and stops, or top quarks and charged Higgs bosons.^] Since SM and 
SUSY contributions appear at the same order of perturbation theory, the measurement of 
the inclusive decay B — > X s ~f given by the CLEO and Belle Collaborations |58[ [as well 



as by the ALEPH collaboration [p9| , albeit with a larger error] is a very powerful tool for 
constraining the parameter space. 

In our analysis we use the value quoted by the Particle Data Group [|15|] : 

BR(6 -> 37) = (3.37 ± 0.37 ± 0.34 ± 0.24+°;^ ± 0.38) ■ 10" 4 . (22) 



16 We neglect contributions from loops involving gluinos or neutralinos, which are much smaller than the 
chargino loop contribution in mSUGRA type models p7j. 
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The first three errors are, respectively, the statistical error, the systematical error and the 
estimated error on the model used to extract information about quark decays from B meson 
decays. The fourth error is due to the extrapolation from the data, where the photon energy 
is cut off at 2.1 GeV, to the full range of possible photon energies which contribute to the 



integrated partial width [60]. The fifth error is an estimate of the theory uncertainty. 

Recently, the next-to-leading order QCD corrections to the decay rate have been cal- 
culated in the MSSM [BTL B3]. In the present analysis, we will use the most up-to-date 
prediction of the b — > S7 decay rate [^], where all known perturbative and non-perturbative 
effects are included.^ This includes all the possibly large contributions which can occur at 
NLO, such as terms oc tan/3, and/or terms containing logarithms of tjisvsy /Mw- Besides 
the fermion and gauge boson masses and the gauge couplings discussed previously, we will 
use the values of the SM input parameters required for the calculation of the rate given in 



Ref. |63 |, except for the cut-off on the photon energy, E^ > (l — 5)mf ) /2 in the bremsstrahlung 



process b — > s'yg, which we fix to 5 = 0.9 as in Ref. |J2 



To be conservative we add all the experimental and theoretical uncertainties in eq. (p2|) 
linearly; note that most of these errors are systematic or theory errors, which do not obey 
Gauss statistics. We thus allow the predicted b — > S7 decay branching ratio to vary within 
the rangeQ 

2.0 x 10~ 4 < BR(6 -> 57) < 5.0 x 10 -4 . (23) 

v) The contribution to the muon g — 2: 

Very recently, the Muon (g — 2) Collaboration in Brookhaven has reported a new measure- 
ment of the anomalous moment of the muon [[1(J : 



(^-2) = a c ^ p = 11659 202(14)(6)10- 10 , (24) 

where the first uncertainty is statistical and the second systematical. The value predicted in 
the SM, including the QED, QCD and electroweak corrections is: af* = 11 659 160 (7) 10~ 10 



|B"5H [see however Ref. [|I8| [19| for the size of the hadronic uncertainties]. This value differs 
from the new average value by 2.6a [pl|. This leads to a 1o range for interpreting the 
discrepancy as a New Physics contribution of: 



11 • 10~ iU < a M (New Physics) < 75 • 10~ iU . (25) 



17 We thank the authors of Ref. |62], in particular Paolo Gambino, for providing us with their FORTRAN 
code and for their help in interfacing it with our code for the MSSM spectrum. 

18 Our conservative approach comfortably accommodates the uncertainty in BR(& — > s-f) related to the 
proper definition of heavy quark masses recently discussed in Ref. |34| . Another argument for a conservative 
interpretation of the bound on BR(& — > 57) is that small modifications of the GUT scale boundary condition 
[specifically, small non-vanishing values for A 2 ^ or A^ 2 ] could greatly alter the prediction for this branching 
ratio, without affecting any of the other quantities discussed in this paper. 
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The contribution of SUSY particles to (g M — 2) is mainly due to neutralino-smuon and 
chargino-sneutrino loops [if no flavor violation is present]. In mSUGRA-type models, the 
contribution of chargino-sneutrino loops usually dominates; it is given by |66||: 



Ao susy = arri^y 



^m^9k9k ( _ 4 1 o 2 

47T 



fc=l,2 

m 



-^77^ vg 1- o^+ q4+ g mx fc 

m~(l — Xfc)' 5 V 3 3 3 



+ — 



fcfc) +(9k 



R^ 2 



3 



3mi(1 _ Xfe)4 v l + ^-3^ + -4 + 3x fc lnx fc 



(26) 



with Xk = m 2 ±/m 2 . The chargino-lepton-slepton couplings are determined by the matrices 

U, V which diagonalize the xf 1 mass matrix: g\ = Uk2/(V^M w sw cos/5) and gj? = —Vki/sw 
|67]. If the SUSY particles are relatively heavy, the contribution of chargino-sneutrino loops 



x 



can be approximated by [rh is the mass of the heaviest particle]: Aa* " ~ 1.3 • 10~ 5 
(tan /3 /in 2 ), to be compared with the contribution of neutralino-smuon loops, Aa* & ~ 
1.1 • 10 -6 x (tan j3 /in 2 ), which is an order of magnitude smaller. 

Thus the SUSY contribution to (g — 2) is large for large values of tan /3 and small values of 
m and mi/2, e.g. reaching the level of 40 • 10~ 10 for tan/3 ~ 50 and rh ~ 400 GeV. Note also 
that the sign of the SUSY contribution is equal to the sign of fi, a^ USY oc (a/7i)fj,M 2 tan (3/rh 4 . 
If the 2.6a deviation of the measured (g^ — 2) from the SM prediction is to be explained by 
SUSY, the sign of \i therefore has to be positive, // > 0. Intriguingly, this sign of /x is usually 
also favored by the constraint ( |23"1) on the BR(6 — > sj) |57|, |62 . 



3.2 Cosmological constraints from the LSP relic density 

In this section, we will discuss the contribution of Xi particles to the overall matter density of 
the Universe, following the standard treatment |68| , with the modifications outlined in |JD| . 
This treatment is based on the following assumptions: i) Xi should be effectively stable, 
i.e. its lifetime should be long compared to the age of the Universe. This assumption is 
natural in models with conserved R-parity if the LSP resides in the visible sector, ii) The 
temperature of the Universe after the last period of entropy production must exceed ~ 10% 
of m^o. This assumption is also quite natural in the framework of inflationary models, given 
that analyses of structure formation determine the scale of inflation to be ~ 10 13 GeV [68 1 . 



However, one should bear in mind that one can construct models where one of these 
assumptions is violated, without changing the collider phenomenology of the (mSUGRA) 
model we are considering. The cosmological constraints we are about to derive therefore 
have a different status than the constraints that follow from "New Physics" searches at 
colliders. It is nevertheless interesting to delineate the regions of parameter space where the 
model can provide the approximately correct amount of Dark Matter in the Universe under 
the stated simple assumptions. 
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If these assumptions hold Xi decouples from the thermal bath of SM particles at an 
inverse scaled temperature xp = m^o/Tp which is given by |68| 

x F = ■ ^-. (27) 

y/9*X F 

Here v is the relative LSP velocity in their center-of-mass frame, cr ann denotes the LSP 
annihilation cross section into SM particles, (...) denotes thermal averaging, Mp\ = 2.4 • 
10 18 GeV is the (reduced) Planck mass, g* is the number of relativistic degrees of freedom 
(typically, g* ~ 80 at Tp), and c is a numerical constant, which we take to be 0.5. One 
typically finds xp ~ 20 to 25. Today's LSP density in units of the critical density is then 
given by |6^ 



., 2.13 • 10 8 /GeV 



O,o/i 2 = — — ' - , (28) 

where the annihilation integral J is defined as 

j (XF)= r^m dx . m 

Jxjp x 

The quantity h in eq. ( p8|) is the present Hubble constant in units of 100 km/(sec-Mpc). 
Eqs. (p7|)(|29|) describe an approximate solution of the Boltzmann equation which has been 



shown to describe the exact numerical solution very accurately for all known scenarios. 

Since neutralinos decouple at a temperature T F <C "rriyO, in most cases it is sufficient 
to use an expansion of the LSP annihilation cross section in powers of the relative velocity 
between the LSPs: 



w a 



w(x\x\ -» SM particles) = a + bv 2 + 0(v A ). (30) 



The entire dependence on the parameters of the model is then contained in the coefficients 
a and b, which essentially describe the LSP annihilation cross section from an initial S- 
and P-wave, respectively. The computation of the thermal average over the annihilation 
cross section, and of the annihilation integral eq. (|29"D , is then trivial, allowing an almost 
completely analytical calculation of fl^oh 2 [eq. (|27|) still has to be solved iteratively, but this 
iteration converges very quickly]. Expressions for the a and b terms for all possible two- 
body final states are collected in ||70| . In these expressions we use running quark masses 



at scale Q = rn^o, in order to absorb leading QCD corrections. Of course, we also use the 
loop-corrected Higgs boson spectrum. 

In generic scenarios the expansion eq. (|30|) reproduces exact results to ~ 10% accuracy, 
which is quite sufficient for our purpose. However, it has been known for some time |jj| that 
this expansion fails in some exceptional cases, all of which can be realized in some part of 
mSUGRA parameter space: 
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- If the mass splitting between the LSP and the next-to-lightest superparticle NLSP 
is less than a few times Tp, co-annihilation processes involving one LSP and one NLSP, 
or two NLSPs, can be important. This can usually still be treated using the formalism of 
eqs. (PTD— fl30"|) if the x?X? annihilation cross section is replaced by a weighted sum of terms 
describing the annihilation of two superparticles into SM particles p9fl . We include the O{v ) 



"a— "terms for co-annihilation with a charged slepton (e^, JIr or f\) [[H|] or scalar top t\ 

[ |7^] . In these cases co-annihilation can reduce the relic density by one order of magnitude 

7% or more J7^] . If xl is higgsino-like, we include co-annihilation with both x\ an d xt [III 



assuming SU(2) invariance to estimate co-annihilation cross sections for final states with 
two massive gauge bosons V from o~(xiXi — ¥ VV). Since LEP searches imply m^o > M w 
for higgsino-like LSP, so that <j{x\x\ ~* W + W~) is large, co-annihilation in this case 
"only" reduces the relic density by a factor < 3. In our numerical scans co-annihilation 
with fi is important near the upper bound on mi/2 for fixed tuq, which comes from the 
requirement vfif x > ^0. Higgsino co-annihilation can be important in the "focus point" 
region m^ 3> rny 2 - Finally, co-annihilation with t\ can be important in some scenarios with 
large \Aq\. 

- The expansion eq. (|30D breaks down near the threshold for the production of heavy 
particles, where the cross section depends very sensitively on the center-of-mass energy 
yfs. In particular, due to the non-vanishing kinetic energy of the neutralinos annihilation 
into final states with mass exceeding twice the LSP mass ("sub-threshold annihilation") 
is possible. We use the approximate treatment of Ref. [|69| for annihilation into W + W~ 
and hh pairs, which can be important for relatively light higgsino-like and mixed LSPs, 
respectively.^ The integral (^) then has to be computed numerically. In mSUGRA, these 
effects can be of some importance for ml 3> rn\i 2 - Note, however, that the bound m^o > Mw 
for higgsino-like LSP implies that sub-threshold annihilation into W + W~ pairs is no longer 
an issue. 

- The expansion eq. ([30]) also fails near s— channel poles, where the cross section again 
varies rapidly with y/s. In the MSSM this happens if twice the LSP mass is near Mz, 
or near the mass of one of the neutral Higgs bosons. In this case we follow the general 



procedure described in Ref. [pSH , using the numerical treatment outlined in Ref. |7^| . In 
mSUGRA, the Z and h poles occur at low mm; most of the Z-boson pole region is now 
excluded by chargino searches at LEP. For small and moderate tan/3 the A and H poles 
are not accessible. However, for large tan/3, i.e. large bottom Yukawa coupling, both soft 
breaking squared Higgs masses are reduced from their GUT-scale input values, so that Ma 



can become quite small [75]. Note that the pseudoscalar A boson pole is accessible from an 



S-wave initial state, while the CP-even H boson pole is only accessible from the P-wave; 



19 Since for m^o ~ Mz annihilation into W + W pairs is already open, a careful treatment of ZZ final 
states is less important; note also that far above threshold, the cross section for W + W~ pairs is about two 
times higher than that for Z pairs. Similarly, we did not encounter situations where a careful treatment of 
final states containing at least one heavy Higgs boson H 1 A or H^ is important. 
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the contribution from if —exchange is therefore suppressed by a factor of v 2 . We will see 
below that the effect of resonant A— exchange can be quite dramatic. 

Recent evidence suggests [|76| that Q^o ~ 0.3 with h 2 ~ 0.5. We define 

0.1 < Q^oh 2 < 0.3 (31) 

as the "cosmologically favored" region. To be conservative we also discuss the region of 
parameter space where 

0.025 < Vt^h 2 < 0.5; (32) 

the lower bound comes from the requirement that Xi should at least form galactic Dark 
Matter, and the upper bound is a very conservative interpretation of the lower bound on the 
age of the Universe. 

3.3 Results 

Using the theoretical, experimental and cosmological constraints discussed in the previous 
sections, we perform a full scan of the (m 1( / 2 , mo) plane0 for given values of the parameters 
tan/3 and Aq, fixing the higgsino parameter \x to be positive to comply with the (g^ — 2) 
constraint. The results are shown in Figs. 1-7, which show the regions in the (mi/2,mo) 
plane excluded or favored by the various constraints discussed above. In Figs. 1-4, the SM 
input parameters and the EWSB scale are as discussed in section 2.2. In Figs. 5 and 6 we 
show the effects of the uncertainties of the top and bottom quark masses, and of the residual 
scale uncertainty, respectively. Fig. 7 shows the effects of the radiative corrections to the 
fermion and SUSY particle masses. 

Let us first exhibit the effects of the individual constraints on the (mi/2, Wo) parameter 
space for tan/3 = 40 and Aq = 0; Fig. 1. The most stringent theoretical constraint, shown 
in the top-left frame of Fig. 1, is the requirement of proper symmetry breaking. In the 
small green area the pseudoscalar Higgs boson mass takes tachyonic values. The region with 
tachyonic sfermion masses is indicated in dark blue. In the yellow area the iteration to 
determine |/x| does not converge to a value fi 2 > 0. The latter constraint plays an important 
role and excludes, depending on the value of tan/3 [and m t , rrib, Mewsb as will be discussed 
later], many scenarios with m ^> mm. The requirement that the LSP is indeed the lightest 
neutralino rules out the region (in light blue) of small values of mo where the less massive 
f\ slepton is lighter than Xv Turning to the experimental constraints on SUSY particle 
masses, the requirement that the lightest charginos are heavier than ~ 104 GeV (brown 

20 For each value of tan/3 and Aq, we vary niQ from 10 to 2500 GeV with a grid of 10 GeV and TO1/2 
from 5 to 1250 GeV with a grid of 5 GeV. This makes 62.500 points for each choice of tan/3 and Aq. The 
maximal value thq — 2.5 TeV corresponds to first and second generation sfermion masses mr > 2.5 GeV, 
while ffij/j = 1.25 TeV leads to gluino masses around 2.75 TeV and squark masses above 2.3 to 2.5 TeV. 
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area) extends the region of no EWSB^], while the requirement of heavy enough sleptons, 
nil ~ 100 GeV, (dark area) slightly extends the region where sfermions are tachyonic. 

The constraint from the measurement of the b — *■ sj branching ratio excludes only a 
small additional part of the parameter space with low m and mu2 values (medium grey 
area) leading to light charginos and top squarks [the constraint would have been stronger for 
jjl < 0]. For our choice \i > mSUGRA generally predicts this branching ratio to be smaller 
than in the SM. However, we will see later that exceptions to this rule can occur for large 
and negative values of Aq. The lightest Higgs boson mass constraint Mh > 113 GeV (in 
both the medium and the light grey areas) is only effective if m < 1 TeV and mi/2 ^ 300 
GeV since we are in a large tan/3 scenario where Mh can easily be sufficiently large. Note 
that for the values of tan/3 and A used here, there are no points, not already ruled out 
by the constraints on EWSB and the SUSY particle mass bounds, which are excluded by 
the Ap constraint, since the splitting between the top squarks remains moderate. The CCB 
constraint, which is somewhat related, is also not effective in this case, because A t (M EW g B ) 
remains moderate compared to the masses of the stop eigenstates. 

Let us now come to the positive indications for SUSY which are shown in colors in 
the three other panels of Fig. 1. The "Obelix menhir" (dark red) in the top-right frame 
corresponds to the region where the lightest h boson mass lies in the range 113 GeV < Mh < 
117 GeV. It extends from mi/ 2 values below 550 GeV for mo ^ 1 TeV down to rai/2 ~ 300- 
200 GeV for larger mo, with upper contour at mo = 1.7 TeV. For larger values of m and 
rrii/2 the top squarks are very heavy and push Mh beyond the 117 GeV limit. One notices 
that the Mh constraint is satisfied in a large region of mSUGRA parameter space, since a 
variation of the h boson mass of a few GeV leads to a variation of mi/2 of several hundred 
GeV. This is due to the logarithmic dependence on the stop masses, which in turn are mainly 
driven by rai/2- Since the theoretical error on Mh [from higher order loop corrections, as 
well as from a shift of M t by a few GeV within the experimental error, etc.] are expected to 
be of the order of a few GeV, we display for illustration the effect of including an additional 
uncertainty of ±2 GeV on the h boson mass which varies then in range 111 GeV < Mh < 119 
GeV (light red area). As can be seen, the impact is very large and values mo ~ 2 TeV and 
m i/2 < 0.8 TeV could be reached [this point has also been realized recently, see Ref. [[ft| for 
instance]. Some caution is therefore needed when analyzing the consequences of a "115 GeV 
h boson" for the allowed mSUGRA parameter space. 

In the blue regions of the bottom-left frame of Fig. 1 the contribution of SUSY particles 

21 For small values of mo the right-hand side of this boundary does not depend on mo; in this region, Xi 
is wino-like and its mass is approximately given by m~± ~ Mi ~ 0.8mi/2. For larger values of mo, one 

enters the "focus point" J77J region where X\ is a mixture of higgsino and gaugino states; for even larger 
values of mo, fi becomes smaller and the chargino is higgsino-like with m~± ~ \/j,\, until one reaches the "no 
EWSB" region where no consistent value of /j, is obtained. Note that close to this boundary, the conservative 
experimental bound m~± > 84 GeV would have been more appropriate, but the strip where the constraint 

would have been different is very small. 
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Figure 1: Constraints on the (mi/2,Tno) mSUGRA plane. Top-Left: individual constraints 
from non-convergent ji (yellow region), tachyonic Ma (green), tachyonic sfermions (blue), 
light sfermions (dark), light charginos (brown), Xi non-LSP (light blue), BR(6 — > s'j) 
(medium grey) and light h boson (light and medium grey). The three other plots are for the 
ler (dark colors) and 2a (light colors) "evidence" for, the Higgs boson (but with larger error 
bars, Top-Right), (g^ — 2) (Bottom-Left) and the Dark Matter (Bottom-Right). 
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to the anomalous magnetic moment of the muon is within two standard deviations (light 
blue) and one standard deviation(dark blue) from the central value of the measurement made 
by the Brookhaven experiment, eq. (|25|). The 2<r area extends from values mo < 1.2 TeV for 
small mi/2 to the boundary where the neutralino x\ is not the LSP for large values, rrii/2 ~ 1 
TeV, except for a little corner with m^ ~ m ~ 250 GeV, where the SUSY contribution 
exceeds the 2<r upper bound. In this area, charginos and smuons have relatively small masses 
and can give too large a contribution to (g^ — 2). The la constraint is significantly more 
severe: a large amount of the upper part of the la area [where charginos/smuons are too 
heavy to contribute] and a smaller area of the lower part [where charginos and smuons are 
too light and generate too large a contribution] are cut away. 

Finally, the light green bands in the bottom-right frame correspond to the regions where 
the LSP neutralino cosmological relic density is in the required range, 0.025 < Q^ah 2 < 0.5. 
The narrow band slightly above the non Xi LSP boundary is the region where fi is almost 
degenerate with the neutralino x\i an d ^1X1 as wen as ^iTi co-annihilation is efficient 
enough to reduce the relic density. The region below m\ji < 400 GeV and small m is the 
bino-like LSP region, where both the LSP and the f\ are light enough for the annihilation 
XiX? - > t + t~ cross section, through i-channel f\ exchange, to be sizeable. The area near 
the EWSB boundary is again the focus-point region where the Xi has a significant higgsino 
component, i.e. has sizeable couplings to massive gauge bosons, making the annihilation 
into WW, ZZ efficient. Very close to the "no EWSB" region the xi is higgsino-like and is 
therefore almost degenerate with the lightest chargino and next -to-lightest neutralino, in 
which case higgsino co-annihilation takes place. The relic density can then even fall below its 
lower bound. Note that we did not yet reach the regime where the s-channel pseudoscalar 
Higgs boson poles play an important role, although we have a relatively large value of 
tan/3. We also show the regions where the cosmological relic density is more constrained, 
0.1 < Q^oh 2 < 0.3 (dark green). The areas become narrower, in particular the bino-like 
Xi and the focus point regions, but there is no qualitative change from what was discussed 
above. 

Fig. 2 summarizes the situation when all the constraints as well as the 2a "evidence" 
from the LEP2 lightest h boson, the (g^ — 2) deviation and the neutralino cosmological 
relic density are superimposed. [Note that some of these colored areas overlap as can be 
inferred from Figs. 1; we refrained from allocating different colors for these common regions, 
since one can deduce them by continuing the boundaries up to the dark region which is 
the intersection of all three areas]. As can be seen, the region excluded by theoretical and 
experimental considerations is still relatively modest. There are large areas of the parameter 
space where one can accommodate a ~ 115 GeV h boson and a SUSY explanation of the 
(g M — 2) deviation. On the other hand, the area where the neutralino LSP is a good Dark 
Matter candidate is fairly small for this value of tan/3. The areas (in black) where all of 
the three requirements are met are rather tiny and include only a part of the region with 
a light bino LSP neutralino and a very small part of the focus point region, the remaining 
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Figure 2: Constraints on the (mi/2,mo) mSUGRA plane for tan/3 = 40, Aq = and 
sign(yu) > 0. The grey areas are those excluded by the requirement of EWSB breaking and 
limits on SUSY particle masses (darker grey), BR(6 — > sj) (medium grey) and Mh > 113 
GeV (light and medium grey). The colors are for the "evidence" for the Higgs boson (red), 
the (g^ — 2) (blue) and Dark Matter (green). 
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pieces being removed by the (g^ — 2) or Higgs boson constraints. Using the more constrained 
scenario with a narrower range for Q^oh 2 and la errors for (g M — 2), would have collapsed 
the overlap region to a narrow strip in the Xi ~ Ti co-annihilation region, with 420 GeV 
< mi/2 < 520 GeV. The lower bound on m 1 / 2 then comes from the upper bound on (g^ — 2), 
and the upper bound on m-i/2 results from the upper bound on Mh- 

In Figs. 3 and 4 we show the same (mo, TO1/2) plane for smaller and larger values of tan/?, 
respectively, choosing A = or — 1 TeV. The most striking features are as follows: 

(i) The region where EWSB does not take place decreases with decreasing tan/3. In 
fact, for tan/3 ~ 5 there is no EWSB problem, and the lower limit on m 1 / 2 comes from the 
requirement m~± ~ 0.75m 1 / 2 > 104 GeV. In contrast, for tan/5 > 50 the EWSB constraints 

Xi ' 

rules out a larger part of the parameter space. For very large values, tan/5 > 60, the 
requirement of EWSB becomes extremely strong and most of the parameter space is ruled 
out. The "no EWSB" region is controlled in part by the top and 5-quark Yukawa couplings 
since the RG evolution of Mh u and Mn d is very sensitive to modifications of At and At, 
respectively. Note that EWSB with tan /3 > 1 requires m 2 H < m 2 H at the weak scale, 
which in turn requires At < \ t , leading to an upper bound on tan/?. The SUSY radiative 
corrections to mj, play an important role here. For positive values of /i, the radiatively 
corrected 6-mass, and thus \b(Mz), is smaller and the problem is postponed to larger tan/?. 
For /i < 0, the situation worsens and one reaches the no-EWSB regime for smaller values 
of tan/3. A resummation of these corrections, as performed in our analysis ||38|| , is therefore 



crucial for an accurate treatment of scenarios with very large tan (5. The no-EWSB region 
is also smaller for sizable (and negative) values of the mixing parameter A which enters 
the RGEs of the Higgs boson masses, the stop/sbottom masses and the SUSY radiative 
corrections to mt, and m t . 

(ii) The region where the Xi is n °t the LSP decreases with smaller tan/3, since the mixing 
in the f sector is proportional to /ztan/3, so that the splitting between the two stau masses 
vn,f x and mf 2 becomes smaller. The net effect is that f\ becomes heavier, for given values 
of m and m^. On the other hand, taking A = — 1 TeV increases the value of \fi\ as 
determined by EWSB, due to the same RGE effect that decreases the "no EWSB" region. 
This increases f mixing and thus reduces m fl , thereby extending the region where Xi is n °t 
the LSP. 

(iii) The b — > 57 constraint is more restrictive for larger tan/? values, although it is still 
not very constraining in the region of parameter space favored by the "115 GeV Higgs boson" 
and the deviation in (g^ — 2). It also becomes stronger^] when A is sizable (and negative), 
since this reduces m^. For example, in the bottom-left panel in Fig. 3 the b —>■ 37 constraint 
excludes the entire region capable of explaining the (g^ — 2) anomaly, so that no overlap 

22 Note, however, the narrow blue strips at small too and small mi/2 in the left panels in Figs. 4. Here 
b — ► S7 decays are "accidentally" suppressed by cancellations between various new contributions. In the 
region below this strip the predicted branching ratio falls above the upper bound of eq.(22), while in the 
excluded region above the strip the predicted branching ratio is too small. 
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Figure 3: Constraints on the (m 1 / 2 ,'^o) mSUGRA plane for values of tan/3 < 40, A = or 
-1 TeV and sign(/z) > 0. The grey areas are those excluded by the requirement of EWSB 
breaking and limits on SUSY particle masses (darker grey), BR(6 — > 57) (medium grey) and 
Mh > 113 GeV (light and medium grey). The colors are for the "evidence" for the Higgs 
boson (red), (g^ — 2) (blue) and Dark Matter (green). 
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Figure 4a: Constraints on the (mi/2, mo) mSUGRA plane for values of tan/3 > 40, A = 
or -1 TeV and sign(yu) > 0. The grey areas are those excluded by the requirement of EWSB 
breaking and limits on SUSY particle masses (darker grey), BR(6 — > sj) (medium grey) and 
Mh > 113 GeV (light and dark grey). The colors are for the "evidence" for the Higgs boson 
(red), and the (g^ — 2) (blue) and Dark Matter (green). 
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Figure 4b: As Fig. 4a, except that we have used the la range for (g^ — 2) and the tighter 
constraint (31) on Q^oh 2 . 
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region exists even at the 2 a level. The Ap [and CCB] constraint also becomes more relevant 
for the small tan/3 and large — A , again due to the possibly large mass splitting between 
the two top squarks. However, the region where Ap is too large is already ruled out by the 
b — > s'-f [and Mh] constraint. 

(iv) For small tan/5 one needs a large mixing in the stop sector [and thus a sizable and 
negative A ] to maximize the h boson mass and to satisfy the bound Mh > 113 GeV. This 
is true, for example, for tan/3 = 5 [top-left frame in Fig. 3]. Even with large stop mixing, 
Mh cannot exceed 117 GeV, so that the whole parameter space is filled by the Mh ~ 115 ± 2 
GeV "evidence" . For increasing tan /5 [which leads to an increase of Mh] , these areas become 
smaller and smaller, until one reaches values tan/5 > 20 where the maximal value of Mh 
reaches a plateau, as is well known. Note that the mixing changes the form of the red area 



since Mh is maximal for \A t \ ~ JQm^m^ [the so-called maximal mixing scenario]; for a 
fixed value of A t , Mh decreases for larger or smaller stop masses which are controlled by m 
and rrii/2. It is curious to note that for tan/5 = 5, Aq < is crucial for allowing a small 
overlap region where all positive "indications" can be explained [at the 2a level, at least], 
whereas for tan/5 = 20 taking Aq = — 1 TeV removes the overlap region that is present for 
Aq = 0; as noted in (ii), the constraint that Xi is the LSP also plays an important role here. 

(v) The area where the contribution of SUSY particles allows an explanation of the 
deviation of the (g^ — 2) from the Brookhaven result is also very sensitive to the value of 
tan/5. For tan/5 ~ 5 and A near zero, the contribution of chargino and smuon loops is 
not sufficient, except in regions excluded by the other constraints. [We stress again that no 
point for \x < can comply with this constraint]. For larger values of tan/5, the (g^ — 2) 
domain becomes larger and for tan/5 > 40, values of 7771/2 or itlq in excess of 1 TeV are still 
compatible [within 2a] with the central Brookhaven result. We also note that (g^ — 2) is 
less sensitive to Aq than the b — > 57 constraint is. In the former case the sensitivity is only 
due to the increased value of |p| required by EWSB, while the latter constraint is directly 
sensitive to A t through the t mass matrix, as explained above. Nevertheless changing A 
from to —1 TeV significantly reduces the values of mo and 777,1/2 required to explain the 
deviation in (g^ — 2), in order to compensate for the increase of |p|. 

(vi) The region where the Xi LSP is required to make the Dark Matter in the universe 
is also very sensitive to the values of tan/5 and A . For small tan/5, the regions with light 
bino-like LSP and the mixed gaugino-higgsino LSP [focus-point] are smaller; the latter is 
even absent for tan/5 ~ 5, or for A = —1 TeV. If tan/5 > 20, this latter choice of A also 
brings the light bino LSP-like region in conflict with the b — > 57 constraint, which becomes 
more severe when A is reduced below 0, as discussed in (iii). [Note that for tan/5 ~ 5, 
there is a region where the required value of fi^oh 2 is attained due to Xi ~ti co-annihilation; 
however, this region is excluded by the h boson mass constraint.] For very large tan/5 
values, the bino-like region becomes much wider, due to reduced f\ mass and larger Higgs 
exchange contributions. Moreover, the mixed or higgsino-like region increases in size, since 
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the focus-point scenario is easier to realize at large tan/3. For tan/3 > 50, there are additional 
domains where Q^oh 2 is in the interesting range, the regions near the pseudoscalar A boson 
or scalar H boson s-channel poles. As discussed previously, for tan/3 3> 1, Ma [and thus also 
Mg\ become smaller, and their Yukawa couplings to b quarks and r leptons [proportional 
to tan/?] are strongly enhanced. The resulting large XiXi ~~ > bb,r + T~ annihilation cross 
sections reduce the relic density to the required level. Note that with our treatment of the 
QCD correction to the bottom Yukawa coupling, these Higgs pole regions open up only for 
values tan/5 > 50 in agreement with the recent analyses performed in Ref. [pL4|| . Moreover, 



the corrections to the physical Higgs masses are of some importance here. The running 
A and H masses are very close to each other in this region of parameter space, but these 
corrections decrease Ma and increase M#0, leading to a mass splitting of 10% or more for 
tan j3 = 60. This is larger than the width of these Higgs bosons, which amounts to typically 
4% of their mass, if —exchange still does not lead to a separate favored region in our scans, 
due to the P— wave suppression, but the large Higgs mass splitting increases the width of 
the cosmo logically favored region. As a result, even the la overlap region becomes quite 
sizable for tan/3 > 50 and Aq = 0, as shown in the right panels in Fig. 4b. 

We now discuss, for the choice tan/3 = 40 and A = 0, the effect of using different top and 
bottom quark mass input values, approximately 1 a and 2<r higher or lower than the central 
experimental values, M t ~ 175 ± 5 GeV and m&(m&) = 4.25 ± 0.125 GeV; see Fig. 5. For 
smaller top quark mass, M t = 170 GeV, the region where EWSB does not occur becomes 
much larger. The M^ constraint becomes also much more severe; one needs significantly 
heavier stops, and hence larger values of mo and/or m^, to obtain a sufficiently large value 
of M/j. The (g^—2) domain remains almost the same, since in the relevant region of parameter 
space the value of |/x| required by EWSB is only slightly reduced by this reduction of the 
top mass. On the other hand, the DM region, and in particular the mixed higgsino-gaugino 
region, becomes wider. For larger M t the trend is reversed: the region where EWSB does 
not occur and the one with large higgsinos-gaugino mixing almost disappear. The Higgs 
mass constraints become less severe, while the domain where 113 GeV < M^ < 117 GeV is 
reduced. 

Since for the features related to EWSB an increase of nib is more or less equivalent to an 
increase of tan/3, our results for mj = 4 GeV and rhb = 4.5 GeV look similar to those where 
the central value nib(nib) ~ 4.25 GeV is kept, but with an input value tan/3 ~ 37 and ~ 44, 
respectively. One sees that in this case, the changes are rather modest. These modest changes 
however conspire to produce a significant reduction of the overlap region when rrib(rrib) is 
reduced from from 4.5 to 4.0 GeV (bottom panels in Fig. 5). The effect of varying nib would 
have been more striking for large tan/3 > 50, where EWSB and the cosmological relic density 
become much more sensitive to the b- Yukawa coupling. For instance, the regime where it is 
difficult to realize EWSB is reached for lower values of tan/?, tan/3 ~ 55, if rrlbijnb) ~ 4.5 



3 Note that A cannot couple to two equal squarks, e.g. ii or b\ pairs, while H does have such couplings. 
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Figure 5: Constraints on the (m 1 / 2 , m ) mSUGRA plane for tan/5 = 40, A = 0, sign(/x) > 
and different values of the pole top quark mass M t and MS bottom quark mass. The notation 
is as in Fig. 2. 
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Figure 6: Constraints on the (m 1( / 2 , m ) mSUGRA plane for tan/3 = 40, A = 0, sign(/x) > 
with different choices of the EWSB scale. The notation is as in Fig. 2 

GeV is used as an input. 

There is a residual scale dependence in our treatment of the MSSM scalar potential, 
since only the one-loop corrections are fully included. A standard RG improvement of the 
effective potential does not completelly remove the scale dependence, due to the presence of 
several a priori unrelated scales. The effect of varying the scale at which EWSB is realized 
[i.e. the scale at which the effective one-loop scalar potential is evaluated and the running 
of the soft SUSY breaking terms is frozen] by a factor of 2 in either direction is displayed 
in Fig. 6. Here again, except for a small change in the shape of the M^ = 115 ± 2 GeV 
domain, the most striking change occurs for the area where the cosmological relic density is 
in the interesting range. Increasing the EWSB scale to M EW sb = 2^/m^m^ (right panel) 
increases the predicted value of \\x\. This reduces the size of the cosmo logically favored 
mixed higgsino-bino region, but leaves the light bino region largely unaffected. On the 
other hand, the choice Mewsb = 0.5^/m^m^ (left panel) leads to a significant reduction 
of the predicted value of |/x|, and a corresponding decrease of the "no EWSB" area as well 
as the cosmologically preferred area where the LSP has a significant higgsino component. 
Moreover, the reduction of |/i| leads to a reduction of the masses of the heavy Higgs bosons. 
This significantly increases the cosmologically preferred region where the LSP is bino-like, 
and hence also increases the overlap region, as can be seen by comparing the left panel in 



24 More sophisticated attempts to further reduce this residual scale dependence have been proposed, see 
e.g. ref. ra. ' 
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Fig. 6 with Fig. 2. In addition a new feature occurs: the opening of the region where the 
s-channel pseudoscalar A boson pole starts to play a role in the LSP annihilation cross 
section. For the present choice of the scale [and tan 0\ this region is still rather tiny, a small 
line parallel to the co-annihilation region. If the scale is reduced to much lower values, e.g. 
to Mewsb = Mt or Mz, this area would have been much more sizeable. However, in this 
case, the requirement of proper electroweak symmetry breaking will exclude large portions 
of the parameter space. 
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Figure 7: Constraints on the (mi/2, m ) mSUGRA plane for tan/5 = 40, A = 0, sign(/x) > 
when certain radiative corrections are not included. In the left panel we have omitted 
corrections to all superparitcle masses, and in the right panel we have in addition neglected 
corrections to heavy fermion masses. The notation is as in Fig. 2. 

Finally, we show in Fig. 7 the effect of the SUSY radiative corrections on the parameter 
space. In the left panel we have switched off all radiative corrections to the masses of SUSY 
particles [neutralinos, charginos, gluinos and squarks] . This has almost no effect on the Higgs 
boson and (g^ — 2) constraints, but the "light" bino-like LSP DM region becomes slightly 
larger. The largest effect is the increase of the region excluded by the b — > 57 constraint: 
the radiative corrections tend to increase the stop and chargino masses, so switching them 
off makes these sparticles lighter, leading to larger contributions to the b — > sj decay width. 



The impact of the SUSY radiative corrections to the top and bottom quark masses [right 
panel] is slightly more important. These corrections tend to increase M t relative to the 
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running mass, so switching them off increases the size of the top Yukawa coupling X t . This 
shrinks the "no EWSB" region, as well as the cosmologically favored region where the LSP is 
higgsino-like. The value of m t that appears in the stop mass matrix, which is the "MSSM" 



running mass in the notation of Ref. [38|, also increases. This shifts the right boundary 
of the area excluded by the b — > s'y constraint to slightly lower values of mi/2- Moreover, 
the region where 113 GeV < Mh < 117 GeV becomes slightly smaller .0 Since the SUSY 
radiative corrections tend also to decrease the bottom quark mass mj, their removal will 
lead to an increase of X b , and hence to a more extended region where the cosmo logical relic 
density is in the interesting range, as is the case for larger mj values [Fig. 5] or larger values 
of tan/3 [Fig. 4]. Finally, we note that switching off all two-loop terms in the RGE would 



dramatically extend the "no EWSB" region, as also pointed out in Ref. |44 



4. Sparticle and Higgs Production in e + e Collisions 

In this section we will discuss the prospects for producing SUSY particles and heavier Higgs 
bosons at future linear e + e~ colliders with center of mass energies between 500 GeV and 1.2 
TeV, in the context of the mSUGRA model. We first list the production processes that we 
will analyze and then discuss the regions of the (m 1 / 2 , mo) parameter space, for various input 
values of tan/3 and Aq, in which these processes are accessible. For completeness, analytical 
expressions for the relevant total production cross sections are given in the Appendix. 

In this exploratory study we will assess the accessibility of certain production modes 
simply through the corresponding total cross section, without performing any background 
studies. However, in most cases the clean experimental environment offered by e + e~ colliders 
should indeed permit discovery of a certain mode, given a sample of a few dozen events. 
Difficulties might arise in some narrow regions of parameter space, which we will point out 
in the following discussion. 

4.1 Production processes 

In our study, we will consider the following production processes: 

Pair production of the lightest charginos : e + e~ — > XiXi"- (33) 

This proceeds through s-channel photon and Z boson exchange as well as t-channel sneutrino 
exchange. For higgsino-like charginos, M 2 ^> |/i|, or heavy sneutrinos, only the s— channel 



25 Recall that our calculation of M/> and BR(6 — » sj) uses M t and rnb(rrib) as input quark masses; these 
are not affected by the SUSY loop corrections. The effect of these corrections is thus entirely through the 
changes of the sparticle spectrum. We also note that switching off the SUSY loop corrections to quark 
masses is not entirely consistent, since we use the routine of Ref. [B2[ for the calculation of BR(6 — > sj) 
which explicitly includes these corrections. We nevertheless feel that Fig. 7 is a reasonable illustration of the 
importance of these corrections. 
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diagrams contribute substantially. Note that the Z boson couples more strongly to the wino 
components of xf than to its higgsino component, and the t-channel v e exchange contributes 
with opposite sign as the s— channel diagrams. The cross section is thus maximal for heavy 
sneutrino and M2 <C |/i|. For small mp e) which corresponds to small mo, the t-channel 
contribution can reduce the cross section significantly. If the chargino is higgsino-like, the 
cross section becomes insensitive to ma, since the ve\\ couplings are zero in this limit. We 
thus see that the production cross section for fixed chargino mass can vary significantly across 
the parameter space, but it is always rather large, so that masses close to the kinematical 
threshold should be accessible. The only exception can occur in the very higgsino-like region, 
where the xf — xi mass difference becomes small. In the most extreme case one may have to 
rely on events with additional hard photon to suppress backgrounds |3(| . However, we found 
that this can happen only in a very narrow strip near the "no EWSB" boundary, where the 
LSP relic density is below the currently favored range. 

The xf decay pattern depends on its mass difference to the LSP, as well as to sleptons and 
charged Higgs bosons. If the lighter chargino is the lightest charged sparticle [generally for 
m o ~ ^1/2] , it mostly decays into Xi plus a real or virtual W boson which in turn decays into 
//' pairs with well-known branching ratios. For smaller ratio mo/mi/2 real or virtual slepton 
exchange contributions become important, in some cases leading to a leptonic branching ratio 
near 100%. There is a very narrow strip in parameter space where Xi ~^ v^ (£ = e, /x, r) 
has a large branching ratio but the charged lepton is very soft. Here one might again have 
to require the existence of an additional hard photon in the event to suppress backgrounds. 
However, in this case charged slepton pair production is also accessible. There are also 
regions where xf decays predominantly into t^u t Xi, either because the fi mass is reduced 
compared to the other slepton masses [this happens at large tan/3 f8l| ], or because L-R 



mixing greatly enhances the xff\v T couplings relative to the corresponding couplings of 
ei ~ eR and fii ~ fin] this latter effect can become important already for moderate values 
of tan/3, if vn\j%j2 < m < fni/2 |82fl - Finally, for very large values of tan/3, charged Higgs 



boson exchange contributions can also play a role, again leading to an enhanced branching 
fraction for the r mode. However, a large or even dominant branching ratio into TV r %\ is not 
expected to significantly degrade the mass reach for xXxi production at e + e~ colliders, which 
should be very close to the kinematical limit [as at LEP]. To the contrary, the measurement 
of xi decay branching ratios might allow one to extract information about (s)particles that 
are too heavy to be pair-produced. 

Production of the lightest neutralinos : e + e~ — > XiX°- (34) 

This process is mediated by s— channel Z— boson exchange and t— and u— channel cl-iGr 
exchanges. In the gaugino limit, |/i| 3> Mi, M2, the Z boson coupling to neutralinos vanishes 
and only the t and w-channel contributions are present. The latter will be suppressed for 
high selectron masses, i.e. for large mo; however, mo ^> mi/2 also generally implies that 
|//| is not so large, so the size of the Z— exchange contribution increases in this region. In 
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the extreme higgsino limit, only the Z boson exchange contribution will survive since the 
eex?2 couplings are (9(10~ 5 ). Except in the extreme higgsino limit, the cross section is 
much smaller than the cross section for chargino pair production; however, as will be shown 
later, the anticipated high luminosity should ensure a detectable signal over most of the 
kinematically accessible parameter space. 

As well known |83], |81[ the x\ branching ratios depend strongly on details of the SUSY 



particle spectrum. For example, the branching ratio into £ + £~Xi can vary between nearly 
100% and almost zero. The former occurs if X2 ~^ £^£i is the only 2-body decay mode 
of x\i while the latter situation is e.g. realized if X2 ~^ ^Xi is dominant. However, at an 
e + e~ collider hadronic X2 decays are as easily detectable as decays into charged leptons. 
The only potentially difficult scenarios are the extreme higgsino region, where the X2 ~ Xi 
mass difference is small [but stays about twice as large as the Xi Xi mass difference] , or 
scenarios where x® decays almost exclusively into the invisible mode vi>x\] however, this 
latter scenario is never realized in mSUGRA, given current experimental constraints. 

We will not discuss the production of pairs of the next-to-lightest neutralinos, e + e~ — > 
X2X2! since in mSUGRA this process leads to the same reach at e + e _ colliders as xt P arr 
production; the approximate equality m^o ~ m~± holds in both the higgsino and the gaugino 
limit. However, the neutralino production cross section is smaller due to the absence of 
the photon exchange channel, and because the Z boson does not couple to neutral SU(2) 
gauginos. Nor will we consider the production of heavier states, e + e~ — > x^X°j with i or 
j > 2 and e + e~ — ■» Xi X.J with i or j > 1, since these channels cannot extend the overall 
discovery reach. Of course, one would eventually like to also study these channels in detail 
in a clean environment. 

Pair production of charged sleptons : e + e _ — * £ + £~ . (35) 

Pairs of SU(2) doublet ["left-'*] and singlet ["right-handed"] selectrons are produced via 
s-channel photon and Z boson exchange and the t-channel exchange of the four neutralinos 
Xi- Since the vector boson couplings to charged slepton current eigenstates are full strength 
gauge couplings and L-R mixing between selectrons is negligible [for our purposes], the 
first two channels always contribute. Since the electron Yukawa coupling is tiny, only the 
exchange of the gaugino-like neutralinos plays a role here. To good approximation the cross 
section is therefore determined by the sizes of the soft breaking gaugino masses Mi and M2. 
In mSUGRA they are related via eq.(l), which leads to Mi ~ M 2 /2 at the weak scale. The 
value of |/i| is not relevant here. Moreover, mixed e"^e^ production is possible through the 
exchange of neutralinos in the t— or u— channel. However, in mSUGRA the overall mass 
reach in selectron pair production will be determined by e^ pair production, since the e^ — e^ 
mass difference can be sizable. 

The production of second and third generation charged sleptons only proceeds through 
s— channel 7 and Z boson exchange. Since the mixing in the f sector can be large, it has 
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to be included in the Zff couplings; in this case we will only consider the production of the 
lighter f states, e + e~ — ► fifi, which offers the largest reach.|^] Note that, with the exception 
of ef e^ production, the cross section for the production of sleptons near threshold is strongly 
suppressed by j3 3 factors, (3 being the cms velocity of the sleptons. Therefore only slepton 
masses up to several GeV below the kinematical limit can be probed. 

The sleptons will mostly decay into their partner leptons and the gaugino-like neutralinos 
and [if accessible] charginos. Slepton pair production at next-generation e + e~ colliders will 
only be accessible if rriQ is not very large, away from the "focus-point" region. The gaugino- 
like neutralinos and charginos will then be the lighter states, since mSUGRA predicts |/i| > 
Mi unless mo ^> mi/2- In particular, for the lighter slepton mass eigenstates the by far 
dominant decay will be if — > £ ± x?- If phase space allows it, the heavier left-handed sleptons 
will decay predominantly into wino-like charginos Xi or neutralinos yQi pl us a neutrino 
or charged lepton, since these decays occur via SU(2) couplings which exceed the U(l)y 
couplings responsible for £ L — >■ £xi decays. In the case of f sleptons, both states might be 
able to decay via the charged current because of L-R mixing. For large tan /3 the T2 decay 
pattern can be rather complicated ||84| . The mass reach for fi pair production is expected to 
be a little lower than that for /} pair production, since the r leptons produced in fi — > rx? 
will themselves decay, which degrades the visible energy. This becomes of some concern in a 
narrow strip of parameter space close to the lower bound on mo, i.e. near the region excluded 
by the requirement m^a < vfif x - On the other hand, f\ pair production generally still gives 
a larger reach in the (mi/2,'^0) plane, since for tan/3 ^> 1, m^ is significantly smaller than 
mjx R . Moreover, the measurement of the polarization of the produced r leptons could yield 
important information about the SUSY model 



Pair production of sneutrinos : e + e — > vv. (36) 

Muon and tau sneutrino pairs are produced only through s-channel Z-boson exchange. 
Electron sneutrinos can also be produced through t-channel diagrams with the exchange of 
the [gaugino-like] charginos, which enhances the cross section significantly if mx/2 is not too 
large. All v pair production cross sections show the familiar j3 3 behavior near threshold. 

In mSUGRA the sneutrinos [as well as the "left-handed" charged sleptons] are usually 
heavier than the SU(2) gauginos. Sneutrinos can thus generally decay through both neutral 
and charged currents into leptons and gauginos, v -^ vx® and v — > I xf- In the narrow 
range where m^o < m^ < m-±,m^o, which can occur for m.1/2 3> m Q if the overall SUSY 
mass scale is rather low so that the negative D— term contribution to m? plays a role, the 
only allowed decay will be the invisible decay into the LSP and a neutrino. However, in 
mSUGRA such scenarios are excluded by the h mass constraint. Usually this invisible decay 



26 The Zfxf\ coupling vanishes for a specific value of the f mixing angle, cos 2 9 T — 2 sin 2 9\y However, 
the photon exchange contribution ensures that the fi pair production cross section remains sizable even in 
this case. 
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is disfavored by the smallness of the U(l)y gauge coupling. 

Pair production of stops and sbottoms : e + e~ — > titx , 6161. (37) 

In e + e~ collisions, squarks can only be produced through s-channel photon and Z boson 
exchange. In our analysis, we will consider only the pair production of the states t± and b\ 
which are in general significantly lighter than the other squark states, due the large Yukawa 
couplings of top and bottom quarks. In both cases L — R mixing can affect the coupling to 
Z bosons, and hence the total cross section p6 , but the latter remains sizable unless these 



processes are kinematically suppressed. 

The mixing in the stop sector and the large top quark mass might generate some com- 
plications in the decays of t\, since the usual two-body squark decays q — > qg, qx® and 
q'Xi are kinematically closed for sizable regions of parameter space. In this case, t\ decays 
are dominated by loop induced two-body decay channels [such as the flavor changing decay 
into a charm quark and the LSP p7|], or by tree-level three-body [bWxi,bH + Xi or Ml] or 



even four-body [bxiff where /, /' are massless fermions] decay modes |5S| . In contrast, 
tree-level two-body decays b\ —* bx® decays are always possible at least for i — 1. As in the 
case of f , the decay patterns of the heavier states 62 and £2 can again be rather complicated 



B4 . 



Pair production of Higgs bosons : e + e — » AH and H + H . (38) 

If the pseudoscalar Higgs boson A is relatively heavy, Ma ?t 200 (130) GeV for moderate 
(large) values of tan/3, we are in the so-called decoupling regime where the running masses 
of the heavier Higgs bosons A, H are almost degenerate^} for even larger Ma, the H ± mass 
also becomes very similar. In that limit, which is realized for most of the allowed mSUGRA 
parameter space, the lighter h boson mass reaches its maximal value, Mh ^130 GeV, and 
h has Standard Model like couplings to fermions and gauge bosons. The pseudoscalar Higgs 
boson has no tree-level couplings to two gauge bosons. In the decoupling limit, the couplings 
of H to two massive vector bosons are also strongly suppressed, as is the ZAh coupling. The 
only significant 2^2 process for the production of the heavier neutral Higgs bosons is 
thus associated HA production through Z boson exchange in the s-channel. The charged 
Higgs particles can be pair-produced photon and Z boson exchange. These cross sections 
are sufficiently large [although once again suppressed by /3 3 factors near threshold] to allow 
for the production of detectable quantities of these particles up to almost the beam energy. 



If tan/3 > JfrT t (fn t ) /frib(fnb) or Ma < 2M t the heavy neutral Higgs bosons will mainly 
decay into bb and t + t~ pairs [with relative abundance ~ 9 : 1]. For smaller values of 
tan/3, decays into ti pairs dominate, if kinematically allowed. The charged Higgs bosons 
will dominantly decay into tb and rv final states. In some cases decays into SUSY particles 

27 We mentioned earlier that even for large Ma, the A and H pole masses can differ by up to 10% if tan/3 
is very large. 
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[charginos, neutralinos and possibly sleptons and stops; only this last mode can compete 
with decays into tt, or with bb at large tan/3] can be possible [ 89f . These decays will be 



more difficult to analyze, but in most cases should be clean enough to be detectable at 
e + e~ colliders. [The main background will probably come from the direct pair production 
of heavier superparticles.] 

Of course, the light scalar h boson can be produced [in association with a Z boson, and/or 
through gauge boson fusion] over the entire mSUGRA parameter space. Detailed studies 
[ |90| , |9l| have shown that detection of this particle is straightforward at an e + e~ collider 



with cm. energy y/s > 250 GeV even with moderate luminosities, / £,dt > 10 fb" 1 . With 
the high luminosities expected at a linear collider, / Cdt ~ 500 fb -1 , detailed studies of the 
profile of this particle can be made already at y/s < 500 GeV: the mass can be measured 
at the permille level, various couplings to fermions and gauge bosons [which determine the 
production cross sections and decay branching ratios] can be measured at the percent level, 
and even the Higgs self-coupling can be measured at the level of ~ 10 — 20%. In the 
MSSM these measurements should reveal some deviations from Standard Model predictions 
for Ma % 500 GeV to 1 TeV. Nevertheless this indirect evidence for the existence of heavier 
Higgs bosons will have to be confirmed through their direct production. 

4.2 Results 

In a first step, we will analyze (s)particle production at an e + e~ collider with y/s = 800 
and an annual luminosity of £ = 500 fb -1 /yr. This is e.g. expected for the second phase 
of the TESLA machine ||92|| . We will consider a given channel to be visible if its total cross 



section exceeds a min = 0.1 fb, which means that a sample of 50 signal events per year will be 
required, with the assumed luminosity [or 100 events with a two years running], to establish 
discovery.^] This number should usually be sufficient for discovery in the clean environment 
offered by e + e~ colliders |9(], |91fl . 

We will also illustrate the potential of the first phase of the e + e~ collider with y/s = 500 
GeV and the same luminosity. We will then discuss the increase of the discovery potential if 
the energy of the collider could be raised to y/s = 1.2 TeV [p5] . This could be accomplished 
either by simply extending the electron and/or positron accelerators, or by increasing the 
RF power if the cavity modules generate higher gradients than originally envisaged [as is the 
case for the most recent TESLA cavities] .0 However, in order to compensate for the 1/s 
drop of most background cross sections with raising cm. energy, we will require a smaller 
minimal cross section than in the previous cases, <r m i n = 0.025 fb. In order to obtain the 
same number of signal events one would then either have to increase the luminosity [as is the 

28 Here, we will only discuss direct production of the new particles. In some areas of the parameter space, 
cascade decays of SUSY particles might allow for the detection of states which are not accessible directly 
since the corresponding cross sections are too small. 

29 We thank Peter Zerwas for a discussion on this point. 
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case for the TESLA machine [[52] where the luminosity is expected to scale with the energy], 
or to extend the running period of the machine. 

Figs. 8-10 show the regions in the usual (mi/2,mo) plane where various superparticles 
and heavy Higgs bosons can be discovered at an e + e~ collider with a cm. energy ^Js = 800 
GeV and a luminosity / Cdt = 500 fb" 1 , for the same values of tan/? and A taken in Figs. 1- 
3. The grey areas are again those excluded by theoretical and experimental constraints. The 
colored regions correspond to the following production processes: e + e~ — > X1X2 (green+red) , 
e + e~ — > xlxi ( re d), e + e~ — > £ + £~ (blue), e + e~ — > vv (purple), e + e~ —> titi (dark blue), 
e + e~ — > bib\ (dark green) and the heavy MSSM Higgs boson production e + e~ — > HA and 
H^H ± (yellow). Note that some of these regions are overlapping. For example, in Fig. 8 the 
chargino and neutralino areas should be extended until the lower boundary where the LSP 
is not the Xij the f region includes the sneutrino and Higgs regions, and the sneutrino area 
includes the Higgs boson region. Moreover, the region where X1X2 production is accessible 
always includes the entire Xi P a i r production region. 

As stated earlier, chargino pair production, e + e~ — > Xi~Xi~; can be probed for mass values 
up to the kinematical limit, m~± ~ 400 GeV for ^/s = 800 GeV. This implies that for wino- 
like charginos, i.e. in the region m ^$1.5 TeV, values mi/ 2 ~ 1.3m^± ~ 500 GeV can be 
probed. For larger values of mo, the chargino becomes a mixture of higgsino and wino, since 
ji is relatively small. This reduces its mass for fixed TO1/2, so that larger values of rai/2 can 
be probed. As long as the chargino is wino-like, the region where Xi P arr production is 
accessible is almost independent of tan j3; however, some parts of this region might be ruled 
out by the experimental constraints, e.g. for tan/5 ~ 5 [top-left frame in Fig. 9] where 
a part has been eaten by the M^ > 113 GeV requirement. As noted in the discussion of 
Fig. 3, choosing Aq = — 1 TeV increases the value of /x required for EWSB, removing the 
higgsino and mixed regions, and the maximal value of mi/2 that can be probed in chargino 
pair production becomes independent of m . 

Searches for the associated production of the lightest and next-to-lightest neutralinos, 
e + e~ — > X1X25 further increase the reach in mi/2 in the bino-like limit. In this region the 
ZX1X2 coupling, which is proportional to the product of higgsino components of Xi an d X25 
is small. However, selectrons are not too heavy, so selectron exchange in the t— channel 
produces a cross section that should be detectable given the anticipated large luminosity. 
The situation is thus quite different from the one at LEP2, where the integrated luminosity 
is three orders of magnitude smaller, in which case the gain from neutralino production 
over the information obtained from chargino pair production is marginal. However, for very 
large values of itlq and A$ = the lighter neutralinos become almost higgsino-like, and 
are therefore almost degenerate with the lightest charginos. The reach in chargino and 
neutralino production then becomes almost the same.Q Finally, for Aq = — 1 TeV the light 

30 We remind the reader that the visible energy released in the pair production of light higgsino-like states 
can be quite small, making the study of these channels difficult in the extreme higgsino region. 
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Figure 8: The (m^, mo) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy i/i = 800 GeV. The grey areas are those excluded 
by theoretical and experimental constraints. The colored regions are those where a given 
cross section is large enough for the particles to be produced: X1X2 (green), Xi~xr ( re d), 
£ + £~ (blue), vv (purple), iiii (dark blue), &1&1 (dark blue) and the heavy MSSM H,A,H ± 
bosons (yellow). Note that some of these regions are overlapping. 
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Figure 9: The (mi/ 2 , m ) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy ^/s = 800 GeV for different values of tan/3 < 40 and 
A with sign(/i) > 0. The notation is as in Fig. 8. 
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Figure 10: The (m 1 / 2 , m ) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy ^Js = 800 GeV for different values of tan/? > 40 and 
Aq with sign(/i) > 0. The notation is as in Fig. 8. 



40 



neutralinos remain gaugino-like even for m > 2 TeV. However, here the selectrons are so 
heavy that the cross section for neutralino pair production becomes too small even at a 
high-luminosity collider. In this region of parameter space even the discovery reach could 
thus still be increased by yet higher luminosities. 

Charged sleptons can be produced for relatively small values of mo- f\ pair production 
remains accessible for values of m from ~ 400 GeV down to the region where the Xi is n °t 
the LSP.Q Sneutrinos are heavier than f±, since they are pure SU(2) doublets. Hence their 
pair production can be probed in a smaller region of the parameter space. Most slepton 
pair production channels are insensitive to tan/?. On the other hand, the region where f\ 
pair production, which offers the largest overall reach in the (mi/2, mo) plane, is accessible 
becomes smaller (larger) with decreasing (increasing) values of tan j3. However, much of this 
changed reach is compensated by the change of the region excluded by the requirement that 
Xi is the LSP. 

If Ao is small or tan/? is large no scalar quarks, not even t\ or b\ squarks, are light enough 
to be produced at ^/s = 800 GeV, given current experimental constraints. At small tan/3 
the Mh > 113 GeV constraint can only be satisfied if squarks are fairly heavy, whereas at 
larger tan/3 the b — > s'y constraint excludes scenarios with light squarks. Only for small 
tan/3 and large (and negative) A does the lightest t\ state become accessible in a narrow 
strip of parameter space with mo ^ 500; see the top-left panel in Fig. 9. 

For tan/3 = 40, Fig. 8, the heavy Higgs bosons H , A and H ± are accessible only in the 
small corner m 1( / 2 ~ 250 GeV, m ~ 200 GeV between the regions ruled out by the Mh 
and the non x\ LSP constraints. This region disappears for tan/3 < 30 [Fig. 9], but can 
be significantly extended for larger tan/3 values [Fig. 10], due to the reduction of Ma with 
increasing tan/3. On the other hand, choosing A$ = —1 TeV rather than increases |/i| as 
determined from EWSB, and hence also the Higgs boson masses. This explains the absence 
of a heavy Higgs region in the top-left panel of Fig. 10. 

By comparing Figs. 1-3 with Figs. 8-10 we see that the entire 2 a overlap region, where 
one has at the same time a light Higgs boson with Mh ~ 115 GeV, a SUSY contribution 
which accounts for the (g^ — 2) deviation and the requirement of a neutralino being a good 
Dark Matter candidate, can be covered by neutralino searches already at an 800 GeV collider. 
Most of this region can also be covered by chargino searches, and much of it can be probed 
in addition by charged slepton (in particular, fi) searches; and sneutrinos are accessible in 
at least part of this region. 

Figs. 11-13 show the SUSY reach for -^i = 500 GeV. The regions where charginos, 
neutralinos and charged (r) sleptons are accessible [with the same requirement for cross 
section and luminosity as before] can essentially be obtained from Figs. 8-10 by simply 
rescaling itiq and m^. The reach is therefore reduced by a factor ~ 5/8 in this case. Of 
course, the experimental constraints remain the same, and thus exclude a relatively larger 

31 Recall that the region very close to the lower bound on m could be difficult to access experimentally. 
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Figure 11: The (mi/2, mo) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy y/s = 500 GeV. The grey areas are those excluded 
by theoretical and experimental constraints. The colored regions are those where some cross 
section is large enough for the particles to be produced: X1X2 (green), x^XT ( rec 0, £ + £~ 
(blue), vv (purple), t\t\ (dark blue), 6161 (dark blue) and the heavy MSSM H, A, H ± bosons 
(yellow). Note that some of these regions are overlapping. 
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Figure 12: The (mi/2, mo) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy y/s = 500 GeV for different values of tan/? < 40 and 
A with sign(/i) > 0. The notation is as in Fig. 11. 
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Figure 13: The (m 1 / 2 , mo) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy y/s = 500 GeV for different values of tan/? > 40 and 
A with sign(/i) > 0. The notation is as in Fig. 13. 
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Figure 14: The (771-1/2, ttio) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy ^/s = 1.2 TeV. The grey areas are those excluded by 
theoretical and experimental constraints. The colored regions are those where some cross 
section is large enough for the particles to be produced: X1X2 (green), x[X\ (red), £ + £~ 
(blue), vv (purple), t\t\ (dark blue), b\b\ (dark blue) and the heavy MSSM H, A, H ± bosons 
(yellow). Note that some of these regions are overlapping. The black line shows the 5<r reach 
contour for sparticles at the LHC in the missing Ey channel with a luminosity / Cdt = 100 
fb _1 ; adapted from Ref. |93 . 
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Figure 15: The (mi/ 2 , m ) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy ^Js = 1.2 TeV for different values of tan/3 < 40 and 
Aq with sign(/x) > 0. The notation is as in Fig. 14. 
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Figure 16: The (m 1 / 2 , m ) mSUGRA plane where SUSY and Higgs particles can be produced 
at an e + e~ collider with a cm. energy yfs = 1.2 TeV for different values of tan/3 > 40 and 
A with sign(/x) > 0. The notation is as in Fig. 14. 
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region of the parameter space that would have been accessible. For example, for tan/3 = 5 
and A = — 1 TeV [top-left frame in Fig. 12] the M^ > 113 GeV requirement has erased a 
sizeable region where charginos or neutralinos can be produced at a 500 GeV e + e~ collider. 
In this case, the stop region has also shrunk dramatically. Moreover, sneutrino and heavy 
Higgs boson production is accessible only in the tiny strip of parameter space at small mo 
and small m-1/2 in the bottom-left frame of Fig. 13 where the new contributions to b — > sj 
decays are "accidentally" small; see the discussion of Fig. 4. 

On the other hand, Figs. 14-16 show that at y^ =1.2 TeV the prospects for SUSY 
particle and MSSM Higgs boson production are of course much better. Again by simple 
rescaling, one reaches values of mi/2 ~ 750 GeV and 1 TeV from chargino pair production 
and mixed neutralino pair production, respectively, while f pair production can be probed 
for m < 700 GeV. The regions where sneutrinos and top squarks can be produced also 
increase, and there is now a little corner where even bottom squarks could be accessible [top- 
left panel in Fig. 15] .0 For large values of tan/?, the Higgs domain increases dramatically 
and for tan/3 > 55 [Fig. 16], the heavy H, A and H ± bosons can be produced in very large 
areas of the parameter space. 

In models with universal soft SUSY breaking gaugino masses the largest value^ of rrii/2 
which can be probed at a 1.2 TeV e + e~ collider through X1X2 production, mi/ 2 — 1 TeV, 
corresponds to squark and gluino masses in the neighborhood of 2 TeV. This is only slightly 
lower than the regions which are expected to be probed at the LHC with a high luminosity, 
where rriq and rrig < 2.5 TeV are accessible p3]| . This is illustrated in Fig. 14, where we 
display, together with the reach of an e + e~ collider with ^/s = 1.2 TeV, the domain that can 
be probed at the LHC with a luminosity of / Cdt = 100 fb _1 in the search for squarks and 
gluinos in the jets plus missing transverse energy channel, shown by the additional line that 
we adapted from the analysis of Ref. [|93|. In addition, there is of course the complementarity 



between the two colliders: while LHC will primarily probe the strongly interacting particles, 
the linear collider will be more concerned with the weakly interacting neutralinos, charginos 
sleptons and Higgs bosons. 

5. Conclusions 

In this paper, we have performed a comprehensive analysis of all the available constraints 
on the minimal Supergravity model. We used the state of the art FORTRAN code SuSpect 
for the calculation of the SUSY and Higgs particle spectrum of the MSSM. We include all 
presently known important effects [two-loop renormalization equations for couplings and 

32 At first sight it is somewhat counter-intuitive that b\ production is accessible only at small tan/3, given 
that mi decreases with increasing tan/3. The reason is that the increase of the region excluded by the 
b — > sj constraint is faster than the increase of the region with small rag . 

33 Note that this value already requires severe fine-tuning to reproduce the experimental values of the W 
and Z boson masses, with typically fi 2 > 10 2 M|. 
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gaugino masses, the complete one-loop effective potential for electroweak symmetry break- 
ing, and the radiative corrections to standard fermion, Higgs boson and sparticle masses]. 
This allowed us to fairly reliably delineate the regions of the mSUGRA parameter space 
which are still allowed by theoretical constraints [from a proper EWSB breaking, neutralino 
LSP, non-tachyonic Higgs and SUSY particles, etc.] and by the present experimental data 
[i.e., the bounds on SUSY particle and Higgs boson masses from negative searches at LEP2 
and the Tevatron, the measurement of the b — > 57 branching ratio, and the precision elec- 
troweak measurements]. We have also indicated the regions of parameter space where one 
obtains a light Higgs boson with a mass Mh ~ 115 GeV, where the SUSY contributions 
can account for the deviation of the experimental value of (g M — 2) from its Standard Model 
prediction, and where the lightest neutralino makes a good Dark Matter candidate. 

When taken at face value, our calculation shows that all constraints [including the positive 
indications] can be satisfied only in fairly narrow regions of mSUGRA parameter space. 
However, experimental errors on the input parameters [in particular, on the mass of the top 
quark; Fig. 5], as well as theoretical uncertainties, are still significant. The latter lead to a 
residual scale dependence of our results [Fig. 6], as well as to a significant uncertainty in the 
prediction of the mass of the lightest Higgs boson. We also repeat our cautionary remarks 
regarding the interpretation of the b — > S7 and, in particular, the Dark Matter constraints, 
which can be circumvented without significantly altering any collider signal, at the cost of a 
somewhat more contrived model of particle physics and/or cosmology. With these caveats 
in mind, the main outcome of our analysis can be summarized as follows: 

- There are large areas of the (Vfti/2,^0) parameter space which are still allowed by 
present experimental constraints. In particular, for large enough values of tan/3, the bound 
on the lightest h boson mass, M h > 113 GeV, does not place too severe constraints. If // > 0, 
which is favored by the (g^ — 2) anomaly, the constraint from the radiative decay b — > s'j is 
not too restrictive [in our conservative interpretation of experimental and theoretical errors] 
even for large values of tan/3. In fact, if A = it is always superseded by the Higgs boson 
mass constraint, but for Aq = — 1 TeV the b — > 57 constraint can be more severe [Figs. 3 
and 4]. Precision electroweak measurements are easily accommodated in the mSUGRA 
framework. 

- For tan/5 > 10 and small values of the trilinear coupling A , the requirement of a Higgs 
boson in the mass range 113 GeV < Mh < 117 GeV favors moderate values of the common 
gaugino mass parameter, mi/2 ^ 500 GeV, leading to relatively light chargino and neutralino 
states, m^± ~ m^o ~ 2m^o < 400 GeV. For large (and negative) values of A , which lead 
to a strong mixing in the stop sector, a Higgs mass in this range can be accommodated in 
large regions of the parameter space even for rather small tan/3(~ 5) values. In this case t\ 
squarks can be rather light, if the soft breaking parameters mo and wiiw are not too high. 
The range of mi/2 favored by the LEP Higgs evidence strongly depends on the exact value 
of M t . We thus eagerly anticipate improved measurements of this important parameter at 
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the current run of the Tevatron collider. 

- The (g M — 2) excess, which can be accommodated in the MSSM only if /i > 0, typically 
gives a stronger upper bound on mo than the requirement Mh = 115 ± 2 GeV. For tan/5 ~ 
40, m and mi/ 2 values below ~ 600 GeV [and slightly above ~ 300 GeV] are needed 
if the experimental value is to be reproduced within la, implying again relatively light 
electroweak gaugino and slepton states. However, the value of this upper bound increases 
roughly proportional to tan/5, so that at tan/3 = 60, m as large as 1.0 (1.6) TeV can be 
accommodated at the 1 (2) a level. 

- For small and moderate tan/5( < 40) the requirement that the density of the lightest 
neutralinos left over from the Big Bang is in the right range to account for the Dark Matter 
density in the Universe is very constraining indeed. In this case most of the region where 
Q^oh 2 is "naturally" in the interesting range is excluded by the Higgs mass constraints, which 
requires SUSY breaking masses above those preferred by Dark Matter calculations. Only 
a small band in the region with a relatively light bino-like neutralino and relatively light 
sleptons survives. In addition, there are "exceptional" regions: a narrow strip in the fix5 
co-annihilation region near the boundary where the fi slepton is the LSP, and a strip in the 
focus point region at large m and small m 1( / 2 values where neutralinos and charginos are 
relatively light and have large higgsino components. Requiring in addition Mh = 115 ± 2 
GeV and a SUSY interpretation for the (g M — 2) anomaly removes most of these "exceptional" 
regions with acceptable relic density. On the other hand, for large values of tan/?( > 50), 
the area of the (m , m^) parameter space favored by cosmology extends significantly due to 
the opening of the pseudoscalar A-boson pole. This allows to fit all the requirements [Mh, 
(g^ — 2) and the DM constraint] in a somewhat larger area of the (mo, m^) parameter space. 
However, the Higgs mixing parameter m\ in the scalar potential needs to be fine-tuned quite 
severely to accommodate very large values of tan /3 |94 . 



- In spite of the strong constraints on the mSUGRA parameter space obtained by taking 
seriously all the positive indications for supersymmetry it is still not possible to give tight 
limits on any one single parameter. We found overlap regions with 5 < tan/? < 60, 0.1 TeV 
< m < 1.5 TeV, 160 GeV < m 1/2 < 550 GeV and for both A = and A = -1 TeV. 
Allowing for a large negative Aq plays an important role in extending the allowed region to 
smaller values of tan j3. On the other hand, for larger tan /3, A = is generally preferred, 
mostly due to the b — > 57 constraint. This implies that the allowed region in the {mi/2,m ) 
plane could be further extended by considering more choices for A , including positive values. 

We have then analyzed the prospects for producing SUSY particles and heavy Higgs 
bosons at high-energy and high-luminosity e + e~ colliders, requiring a sample of 50 events 
per year to establish discovery; this should be sufficient in the clean environment provided 
by e + e~ colliders. At cm. energy yfs ~ 800 GeV, we have shown that charginos, neutralinos 
and sleptons [in particular f sleptons and sneutrinos] are accessible in rather large regions 
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of the parameter space. In particular, already at \fs = 800 GeV associated X1X2 production 
is accessible in the entire overlap region described above. Almost all of this region can 
also be probed through xf pair production, and in much of this region fi pair production 
can also be studied. In some areas of the parameter space, top squarks and even bottom 
squarks can be produced. In the large tan (3 regime, where the present indications for SUSY 
can be accommodated in a larger fraction of the (mi/2,mo) plane, there is a large region 
of parameter space where the heavier MSSM Higgs bosons H, A and H ± are kinematically 
accessible. 

Even for lower cm. energies, \fs ~ 500 GeV, charginos, neutralinos and charged (f) 
sleptons can be produced in a significant region of parameter space not excluded by the 
present constraints. However, discovery of sparticles can then no longer be guaranteed [in 
the framework of mSUGRA] even if all positive indications for SUSY hold up to further 
scrutiny. On the other hand, if the cm. energy of the collider is increased to \fs =1.2 TeV, 
the mSUGRA parameter space where SUSY and Higgs particles are kinematically accessible 
and have sufficiently large cross sections to be detected becomes very wide. The e + e~ collider 
will then have a search potential of SUSY particles that is comparable to the range probed 
at the LHC. This is largely due to the fact that, thanks to the high luminosities expected 
at future e + e~ colliders, the process e + e~ — > X1X2 can probe large values of the parameter 
mi/2' only from kinematical arguments, values rai/2 ~ 1 TeV can be probed at y/s = 1.2 
TeV, corresponding to a gluino mass of the order of 2 TeV. Heavy Higgs particles can be 
searched if their masses are smaller than the beam energy. For large values of tan/3, this 
occurs in a large region of the mSUGRA parameter space. 

One might naively argue that the LHC would still have an edge over a 1.2 TeV e + e~ 
collider, since at the LHC the heaviest sparticles (squarks and gluinos) are produced directly, 
allowing access to lighter sparticles (neutralinos, charginos and sleptons) through cascade 
decays. In contrast, if mi/2 is indeed near 1 TeV [which seems highly unlikely, given the 
original motivation for "weak-scale" SUSY], an e + e~ collider would "only" have access to 
one or two production channels. However, in such a pessimistic scenario the relatively small 
production cross section and the large number of possible decay modes will make it very 
difficult, if not impossible, to study any one decay mode in detail at the LHC. So far at least 



the ability of LHC experiments to perform such studies has only been demonstrated [£3| for 
significantly lower mass scales, i.e. much higher production cross sections. 

On the other hand, precision measurements at an e + e~ collider could reveal a great deal 
about the MSSM spectrum even if only the "lower end" of the sparticle spectrum is accessible 



31]. In particular, thanks to the precise knowledge of the initial e + and e~ beam energies, 
threshold scans allow the measurement of some sparticle masses at the permille level. Making 
use of the ability to vary the beam polarization at will, various couplings appearing in the 
production cross sections of SUSY and Higgs particles can be measured with a high precision. 
Additional couplings can be determined through the careful measurement of decay branching 
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ratios. 

The amount of information on the SUSY model can be maximized by combining results 
from a high-luminosity e + e~ collider with those obtained at the LHC and elsewhere. In the 
near future measurements at relatively low energies [at the B— factories, from searches for 
fi — > ej decays etc.] will yield new information on the flavor structure of the soft breaking 
terms. Any positive signal for large flavor violation would rule out the mSUGRA model 
as commonly defined, while absence of such a signal would strengthen the motivation for 
assuming universal sparticle masses at some high scale. Of course, the crucial test of SUSY 
will have to come from high-energy colliders. A combination of the information on sleptons 
and electroweak gauginos that one can obtain at e + e~ colliders with the information on 
squark and gluino production obtained at the LHC would allow very stringent tests of the 
model. The flavor-conserving parts of the SUSY Lagrangian at the low energy scale can 
then be reconstructed to a considerable extent. Under the assumption of a Grand Desert, 
the underlying structure of the theory at the GUT or Planck scale can then be studied. For 
instance the combination of, on the one hand, the measurement of the electroweak gaugino 
[chargino and neutralino] masses at the linear collider with the measurement of the gluino 
mass at the LHC allows to test gaugino mass unification.^] Similarly, the measurement 
of the slepton and Higgs boson masses at an e + e~ collider together with the squark mass 
determination at the LHC would allow to verify the second mSUGRA assumption that scalar 
soft SUSY breaking mass parameters are also unified at the GUT scale [Q . A high-energy, 
high-luminosity e + e~ linear collider would thus be crucial for fully testing the mSUGRA 
model. 

Of course it is quite possible, perhaps even likely, that mSUGRA in its simplest version 
will already have been excluded before the next e + e~ linear collider commences operation. 
However, most of our conclusions remain valid in a more general supersymmetric context. 
In particular, the SUSY mass scales indicated by the positive evidence for a 115 GeV Higgs 
boson and for a SUSY loop contribution to (g^ — 2) are not very sensitive to details of the 
SUSY model. Our conclusions will then remain qualitatively the same, as long as there 
are no large hierarchies between soft breaking parameters that are assumed to be unified 
in mSUGRA [i.e. as long as all gaugino masses and all scalar masses are similar at some 
high scale; a very large or very small ratio of these two masses can be accommodated even 
in mSUGRA] . In the absence of a compelling model of supersymmetry breaking it becomes 
absolutely essential to collect as many independent pieces of information about the soft 
breaking terms as possible. By studying a highly constrained model we have thus chosen a 
scenario which minimizes the advantages of e + e~ colliders. As argued above, even in this 
context the prospects for SUSY studies at such colliders seem very bright. We are thus 
confident that linear high energy, high luminosity e + e~ colliders will play a crucial role in 



34 A partial test, of the unification of U(l)y and SU(2) gaugino masses, is possible using e + e~ data alone 
pa ]. However, barring a major disaster during construction, data from the LHC will become available first, 
and should thus be included in this test. 
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revealing the secrets of Supersymmetry breaking, assuming that Nature indeed makes use 
of the beautiful idea of weak scale Supersymmetry 
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Appendix: Production cross sections 

In this Appendix, we present for completeness expressions for the production cross sections 
in e + e~ collisions for diagonal and mixed pairs of charginos j9?f , neutralinos |98[], sleptons 
[ P§H and squarks |8B| as well as for the heavy MSSM Higgs bosons ||100| , |5U| , using a unified 
notation. 



1. Neutralino and Chargino production 

The integrated cross section for mixed neutralino XiX? P arr production can be written as 
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with Ajj being the usual phase-space function with the reduced masses: 
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The contribution of the s-channel Z-boson exchange, the t-channel e exchange and the Ze 
interference are expressed in terms of A^ and the final particle energies, e it j = J\j + /i 2 j 
[note that the neutralino masses here are the eigenvalues and the sign must be included]: 
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with Dz{s) = s — Mg + iTzMz- The various couplings are given, in terms of the weak isospin 
and electric charge of the electron and the elements of the matrix Z which diagonalizes the 
neutralino mass matrix [which can be found in Ref. |J7| for instance] by: 
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and 
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For chargino pair production, the cross section can be decomposed into the s-channel 
7, Z-exchange contributions, the i-channel v e contribution and the interference terms: 
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Using the same notation as previously [but now fi 2 = m 2 ±/s and fi 2 = m 2 ±/s] one has for 
the different components: 
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with the additional couplings 0^, O t f expressed in terms of the elements of the matrices U 
and V which diagonalize the chargino mass matrix 
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2. Selectron and Sneutrino production 

The integrated cross section for the pair production of left-handed or right-handed selec- 
trons, which occurs through the s-channel 7 and Z boson exchanges and the t-channel 
exchange of the four neutralinos xf, can be written as 
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The notation is as before with fa is the selectron velocity (3 2 = 1 — 4/x 2 where in this case 
fi 2 = m 2 ./ s. The couplings o» of the selectrons to the Z boson and the couplings \k between 
the electron, the selectrons e, and the neutralinos \°k are gi ven by 
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with the rotated matrix elements: Z' kl = Z kl cw + Z k2 sw and Z' k2 = —Z k i,$w + Z k2 cw The 
kinematical functions F and H, in the case where the two exchanged neutralinos \°i an d X k 
are different, is given by 
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where /i 2 = m|o/s,//? = m 2 - 1 s and the function Fjfc reads: 
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In the case where the exchanged neutralinos are the same [i.e for the squared amplitudes], 
the function Hu k with I = k reduces to: 
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For the production of selectrons of different types, there is no s-channel gauge boson ex- 
change and the cross section simply reads: 
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where in terms of the phase space function defined in eq. ( |A.2| ), 
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For the pair production of the electron sneutrino, the expression of the total cross section 
is similar to the one of left-handed selectrons except that the s-channel photon exchange is 
absent and the couplings are different. It is given by: 
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where the functions H and F are given by the previous equations with /a, replaced by /^ 
and 
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3. Sfermion pair production 

For sleptons of the second and third generation and for squarks, there is only s channel 
gauge boson exchange and the production cross sections, in the absence of sfermion mixing, 
is simply given by: 
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with N c the color factor, N c = 3(1) for squarks (sleptons), and the c^ are as in eq. (|A.11|) for 
a given charge and isospin. However, in the case of third generation sfermions, the mixing 
between the left-handed and right-handed states has to be included. In this case the cross 
section becomes slightly more involved and can be written as: 
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lift. 



na 2 N c 3/ 2 






Ll + Rl 



a%s 2 



1C%rS%r \D Z (S) 



,(A.20) 



with Xij the phase space function eq. ( |A.2| ) and the couplings a^ given by 

a n = Hl cos2 ^/ _ Q/ S w? a 22 = ^at sin 2 ^ - QfSw, a 12 = a 2i = - l[ L sin ^ cos 8j, (A. 21) 

with 6*r the angle of the unitary matrix which turns the left- and right-handed current 
eigenstates into the mass eigenstates: 



/ x = cos 0ff L + sin 9 f -f R , f 2 = - sin 9jf L + cos 9jf R . 



(A.22) 



In the case of squarks, one can include the QCD corrections which can be rather important 
[ |101| , |102| . The standard corrections, with virtual gluon exchange and gluon emission in the 
final state, lead to an increase of the total cross section by ~ 15% far from the kinematical 
threshold, with much bigger corrections closer to threshold. In the case of diagonal pair 



production, they can be included by using the Schwinger interpolation formulae [|101 

r 2 



a(e + e 



QiQi 



a 



Born 



^fe-K^- 6 ) 



(A.23) 



which, up to an error of less than 2%, reproduces the exact results. The corrections for gluino 
exchange are in general smaller and decouple for heavy gluinos. QED threshold corrections 
to slepton pair production can also be of some importance [|103|1 . 



4. Higgs boson production 

The main production mechanisms of neutral Higgs bosons at e + e~ colliders are the Higgs- 
strahlung process and pair production, 



Higgs-strahlung : e + e 

pair production : e + e" 



(Z) -> Z + h/H; 
(Z)^A + h/H; 
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as well as the WW and ZZ fusion processes, 

fusion processes : e + e~ — ► vv {WW} — > vv + h/H ; 

e+e- -► eV(ZZ)^eV + h/H. 

The charged Higgs particle can be pair produced through virtual photon and Z boson ex- 
change, 

charged Higgs : e + e~ — ► (7, Z*) — ► H + H~ 

The production cross sections for the neutral Higgs bosons are suppressed by mixing angle 
factors compared to the SM Higgs production, 

a(e + e- -► Zh) , cr(W -► h) , a(e + e" -► AS") ~ sin 2 (/3 - a); (A.24) 

a(e + e- -> Z#) , <r(FV -> if) , a(e + e - -> A/i) ~ cos 2 (/3 - a), (A.25) 

while the cross section for the charged Higgs particle does not depend on any parameter 
other than Mg±. In the decoupling limit, Ma ^> Mz, the HVV couplings vanish, while the 
hVV couplings approach their SM Higgs values: 

g HVV = cos(/3-a) ^sin4/3M£/2MJ|->0; (A.26) 

g hVV = sm((3-a)^l-0(Ml/M\) -► 1. (A.27) 

Hence, the only relevant mechanisms for the production of the heavy Higgs bosons in this 
limit will be the associated pair production and the pair production of the charged Higgs 
particles. The cross sections, in the decoupling limit and for y/s ^> Mz, are given by [we 

use M H ~ M A ] 

n 2 /\/f 4 
a{e+e-^AH) = ^(^a^l (A.28) 



a(e + e- -> H + H~ 



2G F M w s w 

3tts 



V e VH . (a 2 e + V 2 e )v 



H 



° S W C W ^dUCyySyy 



0h±, (A.29) 



where (3j = (1 — 4M 2 /s) 1 / 2 is the velocity of Higgs bosons, the Z couplings to electrons are 
given by a e — — 1, v e — — 1 + ^s^?, and to the charged Higgs boson by vh = —2 + ^s\y. 

The cross sections for hA and HZ production vanish in the decoupling limit since they 
are proportional to cos 2 (/3 — a). The cross section for the fusion process, e + e~ — > v e v e H, 
is enhanced at high energies since it scales like M^> 2 log s/M\. This mechanism provides 
therefore a useful channel for H production in the mass range of a few hundred GeV below 
the decoupling limit and small values of tan/3, where cos 2 (/9 — a) is not prohibitively small; 
the cross section, though, becomes gradually less important for increasing M# and vanishes 
in the decoupling limit. The cross section for the ZZ fusion process is one order of magnitude 
smaller than that for WW fusion. 
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